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ABSTRACT is 


A discrete-event simulation of a stochastic process, a 
machine-repairman model, has been programmed on the IBM 
Personal Computer. The model consists of three helicop- 
ters, of which two are in service wand one 1s in “coue 
standby, with an option of one or two repairmen. 

The program output 1S a graphics display containing a 
system state-versus-time graph, a table of statistics, and 
animated figures that illustrate the current state of the 
system. The program user can directly observe the dynamics 
of the model as the fixed-increment, simulation clock 
advances. The user has the option of changing the follow- 
ing model parameters: helicopter failure rate, repairman 


Service rate, and the number of repairmen to employ. 
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I. INTRODUCTION 


During recent years, the versatile microcomputer has 
found a niche in many of the professional discliplines 
including education. For difficult and complex subjects, 
educators are exploring the possibilities of the microcom- 
puter aS a teaching aid in addition to its more traditional 
role aS a computing machine. Some universities now use 
microcomputer graphics to illustrate principles that cannot 
be fully explained ina lecture. [Ref. 1] 

A widely used application of computers is to imitate or 
Simulate Stochastic processes. This technique usually 
involves only the numerical evaluation over time of the 


corresponding model with the intent of gaining some under- 


Standing of the process. But uSing graphics as part of the 
output in a Simulation provides the user another 
perspective. For example, observing a working model--such 


aS a birth-death process--can help in grasping the funda- 
mentals of the model and provide some insights into the 
process not easily seen on the classroom chalkboard. 

This paper describes the development of a discrete- 
event, graphic simulation of a specific example of the 
machine-repairman model. The simulation gives the observer 
the impression of real-time passage. A graph on the visual 


display shows the history of the system's transitions from 


state to state, and animated figures illustrate the 
dynamics of the model. The simulation also exhibits a 
table of the theoretical, and the current, estimated values 
of quantities commonly used in evaluating queueing systems. 
Section II of this paper is a discussion of the 
machine-repairman model, the assumptions upon which the 
model is based, and the methods used: for obtaining the 
theoretical and estimated values of the quantities of 
interest. Section III covers the simulation's implementa- 
tion on the IBM Personal Computer (IBM PC). Appendix A is 
a user's guide which contains summary information about the 
machine-repairman model and instructions for operating the 
program. The program flowchart and the program listings 


are found in Appendixes B, C, D, and E. 


II. THE MACHINE-REPATRMAN MODELS 


A. DiSCcvse ren 

The simulation described in this paper actually con- 
tains two machine-repairman models. In both models there 
are three machines, represented- by helicopters. The models 
are distinguishable by the number of repairmen; there is 
one repairman in the first and two repairmen in the second. 

The nominal mission is to keep two helicopters flying 
at all times. The third helicopter is placed in cold 
standby as a backup to the other two. If available, a 
repairman is assigned immediately to a helicopter upon its 
failure. If a repairman is not free, then the helicopter 
joins a queue to await Peon: 

Several assumptions are made to simplify the models. In 
the two-repairmen model, repairmen are assumed to be 
equally competent and to take the same mean time to repair 
a helicopter. Both models assume the repair time to be 
continuous; that is, repair parts are always available, so 
there is no delay in repairing a machine other than the 
time required by the repairman to perform the work. The 
repairmen work independently and do not assist one another. 
Similarly, the models assume the helicopters operate inde- 
pendently and have the same service life with the same mean 


time to failure. By assumption, a helicopter is not 
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subject to failure while in cold standby since it is not 


operating. 
A machine-repairman model typically assumes its 
machines operate continuously until failure. In the case 


of aircraft, this assumption does not reflect reality since 
MP~eEaeeucdl spractice alrcraft do not fly continuously, but 
have instead periods of inactivity between missions of 
limited duration. As already noted, one can reasonably 
assume that an inactive helicopter is not subject to fail- 
ure and that the failure process occurs only when the 
helicopter is active. 

The actual failure rate of a flying helicopter is the 
number of failures per flight hour. The product of the 
actual failure rate times the ratio of the number of flight 
hours performed to the number of calendar hours elapsed 
until failure defines the helicopter's effective failure 
rate. Thus the effective failure rate is the number of 
failures per calendar hour. To simplify the simulation, 
the models assume that the helicopters operate continuously 
until failure with the effective failure rate. 

The selection of helicopters to represent the machines 
in these models allows a simple graphic design in depicting 
their status. An operating machine is seen in the display 
aS a green, flying helicopter; a helicopter that is in cold 


Standby is also colored green but is not flying. A machine 


ia! 


that 1s down and awaiting repair is colored red. When the 
down machine is in repair, it is seen with a repairman. 
Figure l1 is a frequently used, pictorial representation 
(see [Ref. 2], for instance) of the machine-repairman model 
with one repairman. The dashed lines enclose the model's 
repair-queue system which contains a single-seame queue. 
The machines outSide the dashed box are the potential 
customers for the repair queue. Since there are only three 
machines in the model, the longest possible length of the 


queue will be two. 





Figure l: Machine-Repairman Model with One Repairman. 


This model has a second queue, the cold-standby queue, 
Since a helicopter may have to wait before being assigned 
to a mission. The standby queue can contain only one heli- 


copter, and this occurs when there are no machines in the 


repalir-queue system. 


LZ 


The cyclic nature of the model is clearly evident in 
Piemiigqures. A helicopter"s life cycle consists of an oper- 
ating period followed by, after a possible delay in queue, 
a repair period which is in turn followed by, after another 


possible delay in cold standby, another operating period. 


——_ a 
. 
Repalfrman 2 


Figure 2: Machine-Repairman Model with Two Repairmen. 







Standby 
Queue 





Figure 2 depicts a pictorial representation of the 
machine-repairman model with two repairmen [Ref. 3]. In 
this system, helicopters waiting in the repair queue go to 
the first available repairman for service. The maximum 
possible length of the repair queue is only one. As in the 
case of the one-repairman model, the standby queue can 
contain at most one helicopter. 

The machine-repairman model iS a special case of a 
Been—ceath process in which the event of a machine failing 


at time t constitutes a birth. Let T be the operating life 
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of a helicopter until failure and assume that T is exponen- 
tially distributed with a failure rate a. Then the mean 
time to failure (MTTF) of a helicopter is E(T) 23s 
Similarly, the completion of repair on a machine is a 
death, so define S to be the service time of a repairman. 
The service time is also assumed to have an exponential 
distribution with a service rate uw. The mean time for 
repair (MTFR) is then E(S) = 1/u. 

For these machine-repairman models, the state of the 
system, N(t) =n, is defined as the number of helicopters 
down at time t. Since there are only three helicopters in 


the models, then yn = (7-2. 


B. MODEL ESTIMATORS 

In the study of queueing models, there are four quanti- 
ties that are commonly used to evaluate the queueing sys- 
tem's performance. As part of the graphical output, the 
Simulation will provide the user a table of theoretical and 


estimated values of these four quantities: 


L = the average number of machines in the repair 
system, 
L = the average number of machines in the repair 
q queue, 
W = the average amount of time a machine is down, 
W = the average amount of time a machine waits for 
q repair. 
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As shown in Figures 1 and 2, the variable gq is the time 
that a down machine waits in the repair queue. Then 
ue = E(q). It follows that the average total time that a 
helicopter is in the repair system is the sum of the 
average waiting time in queue and the expected service 
time: W = ne Toe uGos a: 


Little's formulas show the relationship of the four 


quantities of interest 


ioe 


eee 
q q 


where \ is defined as the average arrival rate of machines 
entering the repair-queue system. 


There are three additional quantities that will be 


computed and presented by the simulation. They are defined 
as 

eo P{down machine must wait for repair}, 

Ee = P{up machine must go to cold standby}, 

Av = availability = the proportion of time that at 


least one helicopter is not down. 


To calculate all of eight of these values, it is first 
necessary to compute the system's steady-state probabili- 
ties, Pp, = PANo— nei —- 1 0,),2,3. in other words, De. is 


the long-run probability that the system is in state n. The 


details of calculating Por a function of the machine 
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failure rate and the repairman service rate, will not be 
discussed in this paper. The interested reader should 
refer to the bibliography for a list of textbooks on queue- 
ing theory and birth-death processes. In addition to the 
limiting probabilities, {pos n =07 lye the program 
computes the fraction of the time the system is in state n, 
denoted Bae aS equal to the total time in state n divided 
by the total elapsed time. The quantity B, is an estimate 
of Py: 
The value of L, the theoretical average number of 
machines in the repair system, is found by averaging the 
number of down machines over all states: 
3 


t 
ia 
né 


yen = (1)p) 7 (2D (3)p.. 
The estimate of L, denoted L, is computed with the samme 
expression with the exception that the estimate B, is sub- 
SEC ULeds fOr Do: 

The average number of machines in the repair queue Ly 
and its estimate bg depends in part upon the number of 
repairmen in the model. In the one-repairman case, 
machines occupy the repair queue only when the system is in 


State 2 or 3, thus 


L = P5 a3 2P2) and 


cP 
I 


Bp a 2p. - 
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In the two-repalirmen case, a queue forms only when the 


system is in state 3, thus 


The average rate at which helicopters enter the repair 
gqueve, A, 1s dependent upon the number of active machines 
Subject to failure in state n and the proportion of time 
the system is in state n. There are two active helicopters 
in states 0 and 1, one in state 2, and none in state 3. The 
combined failure rate of the active helicopters in state n 


is given by 


= Y n- 2? 
n 
0 a= 
Maen } is found by 
~ 
r Sa ane, = 20D 9 2 Zap, + ADo- 


moi the estimate \, the program divides the number of fail- 


ures that have occurred in time t by t, that is 


X = total number of failures/total elapsed time. 


ley, 


The average waiting time quantities, W and Wo are then 


calculated by uSing Little's formulas 


W = L/A 
W. = le 
q q/ 
To compute the estimate Wa the program keeps a record of 


the length of time each helicopter spends in the repair 


queue and uses the expression 


m 

Wy = (1/m) ) a m = {Opes One 
B= 

where q, is the waiting time in queue of the ith helicopter 

and m is the total number of helicopters that have entered 


the queue by time t. Then W is easily found with 


W=W of Sees 
: cc) 


The estimate PB, of the probability that a down machine 
must wait for repair is the proportion of the time that the 
system spends in a state in which a machine is down with no 


unemployed repairman. The repairman iS not available in 


the one-repairman model when the system is in state 2 or 3: 


In the two-repairmen case, the repairmen are busy when the 


system is in state 3: 
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For ae the only time a machine must wait in cold standby 
Domne weenesystem 1s in state 0; therefore, Eee aioe 

The availibility of at least one helicopter in the 
system is the sum over all states n of the percentage of 
the flyable helicopters in each state multiplied by the 


proportion of time the system is in that state; 


Foe eet ele eee? 7/3) tape el ship (0): 


The estimate of the availabiliy is calculated by substitut- 


ing Doe n = 0,1,2,3, into the above expression. 
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Ill. THE SIMULATION 


A. FH COMPURER 

The microcomputer used for this simulation was the IBM 
Personal Computer. Since the program MACHREPR uses color 
graphics extensively, the computer requires a color display 
monitor and either IBM's Color/Graphics Adapter or a suit- 
able graphics card from another manufacturer. The program 
wlll not run without a graphics adapter. Other minimum, 
hardware requirements for this simulation are 128 kilobytes 
of read-write memory (RAM) and at least one 54-inch, floppy 
disk-~-drive. 

At the beginning of this project, the author chose to 
use IBM's disk operating system PC-DOS with the intent of 
enabling the simulation to run on other microcomputers 
using the more generic version of PC-DOS: MS-DOS by 
Microsoft. It soon became apparent, however, that MACHREPR 
would need to directly access certain parts of the 
computer's memory in order to improve the program's 
performance. Since the program makes direct calls to the 
PC-DOS, there can be no assurance that the simulation will 


run on other "IBM-compatible" machines. 
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B. THE PROGRAMMING LANGUAGES 

MACHREPR is written in Advanced BASIC which is the pro- 
gramming language supplied with PC-DOS. Advanced BASIC was 
selected because of its convenient graphics commands and 
its ready availability. Although Advanced BASIC is. an 
interpreter, program speed was not an essential criterion 
that warranted uSing a compiled language. Advanced BASIC 
was insufficient, however, in two aspects, and so MACHREPR 
uses two subroutines written in 8088 assembly language to 
improve program performance. These two subroutines, named 
SCRNSHFT and RNGEN, will be described later in this 


Section. 


C. PROGRAM DESCRIPTION 

Appendix B contains the milgueneres for the BASIC pro- 
gram MACHREPR. The flowcharts illustrate the simulation's 
general construction and program flow; however, they are 
not a verbatim duplication of the program listing. The 
program listing for MACHREPR and the two source-code list- 
ings for RNGEN and SCRNSHFT are in Appendixes C, D, and &E 
respectively. 

After initializing several variables and arrays, the 
program presents the title screen. The program menu 
appears next with the following options: (1) see the pro- 
gram instructions, (2) change the model's parameters, (3) 


use the default parameters, (4) set the random-number 
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generator seed, or (5) end the program. The instructions 
that are avallable on the screen are a summarized version 
of those found in Appendix A, the user's guide. The model 
parameters that can be changed by the user are the mean 
time to failure, the mean time for repair, and the number 


of repairmen. 


MACHINE - REPATRMAN MODEL 
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Figure 3: An Example of the Graphics Display. 


Figure 3 showS an example of the graphics display 
Screen aS 1t appears during the simulation. The upper half 
of the screen contains the system state-versus-time graph. 
The system state values are on the ordinate of the graph, 
and fifty divisions representing time are on the abscissa. 
The three graphic figures of helicopters in the lower left 


corner symbolize the status of the system. The helicopter 
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seen above the others iS in active service; the two 
machines with broken rotors are down. In this example, the 
helicopter graphics indicate that the system is in state 2, 
where helicopter B is flying, helicopter C is under repair, 
and helicopter A is awaiting repair. The graph shows that 
the current time is 69 and that the system jumped from 
state 1 to 2 at t = 68. 

The table in the lower right corner contains the esti- 


mated and the limiting values of the quantities discussed 


in Section II. The row labels are: 

Ar - the arrival rate of machines into the repair 
queue, 

Wq - the average wait time of a machine in the repair 
queue, 

Lq - the average number of machines in the repair 
queue, 

W - the average amount of time a machine is down, 

L - the average number of machines in the repair 
system, 

Avl - the proportion of time that at least one machine 


1S available. 


Pd - the probability that a down machine must wait for 


repair, 
Ps - the probability that an up machine must go to 


cold standby. 
MACHREPR consists of four major parts: (1) main pro- 
gram, (2) clock module, (3) failure module, and (4) repair 
module. The main program handles the housekeeping details 


such as initializing variables, loading the machine 
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language subroutines, generating the graphics display, and 
performing the tasks selected from the program menu. Within 
the main program are routines for computing the limiting 
values of the tabulated quantities and for printing them on 
the display screen. The main program also calculates the 
initial failure times for the first entry into the clock 
module. 

The clock module determines the next-event time, 
whether the event is a failure or the completion of repair, 
and which helicopter it affects. The simulation clock is a 
fixed-increment, time-advance type. In contrast to a next- 
event, time-advance clock which skips over the intervening 
time between events, the clock module in MACHREPR advances 
the time in equal increments. With this typeof clock >mEm. 
user experiences the illusion of real-time passage. 

During the execution of the simulation, the ellees 
module monitors the keyboard for certain commands issued by 
the user. Pressing the "P" key pauses the program to allow 
the user to study the graphic display. Pressing the "C" 
key causes the program to continue. The "S" key will stop 
the simulation and return the program to the program menu. 
Should the clock reach the preset limit of 9950 time units, 
the program will automatically stop the simulation and 
return to the program menu. 

Every time the clock is advanced one increment, the 


clock module updates the system state-versus-time graph 
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according to the current state and time values. As _ the 
time advances, the plotted lines will eventually reach the 
edge of the graph. When this condition is met, the clock 
module shifts the entire graph left one time unit and draws 
the new time increment. This procedure continues for all 
successive time increments. When the simulation clock has 
reached cnet time for the next event, the program jumps to 
either the failure module or the repair module as 
appropriate. 

The failure module changes the helicopter display 
according to the current state and model parameters. The 
program keeps track of the machines individually, so that 
the graphic image of the particular, failing helicopter 
transforms from a green, flying figure to a red, grounded 
one. If a repairman is free, the graphic figure of a man 
1s drawn next to the aircraft. If no repairman is avail- 
able, then the human figure is omitted, signifying that the 
helicopter is placed in the repair queue. If a repairman 
was asSigned, the failure module calls the random-number 
generator and computes the service time. The sum of the 
service time and the current time is the time that the 
repair will be completed and this is compared to the next- 
event times of the other machines after the program returns 
to the clock module. 

In a similiar manner, the repair module changes the 


graphic figure of the helicopter that has just completed 
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repairs. If there are less than two helicopters flying 
when the repair event occurs, then the repaired machine 
transforms from red colored and disabled to green colored 
and flying. If two helicopters are already flying; eeemem 
the repaired machine goes into the standby queue. The now 
free repairman moves to the helicopter at the head of the 
queue, or if the queue iS empty, he disappears from the 
screen. 

Before the program returns to the clock module, the 
repair module computes the interarrival time to failure of 
the helicopter just repaired. Failure times are not com- 
puted for helicopters placed in the standby queue. If the 
repairman iS assigned to another helicopter, the repair 
module also computes the service time. 

After returning to the clock module, the simulation 
updates the estimates of the tabulated quantities and 
prints them on the display screen. Both the failure and 
repair modules change the state-indicator variable to its 
new value so that the clock module will draw a horizontal 
line one time unit in length at the proper place Oona 
graph. 

SCRNSHFT is the subroutine that shifts the system 
State-versus-time graph one time unit to the left per clock 
cycle. The subroutine uses certaln program transfers 
called interrupts to the PC-DOS Basic Input/Output System 


(BIOS) which contains routines that control the ~vrdee 
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display. SCRNSHFT is not appreciably faster than a BASIC 
routine designed to perform the same task; however, the 
subroutine eliminates an annoying flicker in the display 


that was characteristic of the BASIC version. 


D. PROGRAMMING CONSIDERATIONS 

The program computes the interarrival (failure) and 
service times, assumed to be continuous, exponentially dis- 
marDubked, randem Variables, in the single-precision format. 
The current time, an integer value, is added to the inter- 
arrival time and to the service time to yield the next- 
failure time and the next-repair time respectively. In the 
unlikely event that the two  single-precision, random 
variables match, the simulation will proceed as if the 
failure occurred earlier. 

Whenever the model jumps from one state to another, the 
program computes a next-failure time and a next-repair time 
for the helicopter just affected by the transition. Thus 
MACHREPR always keeps in memory six next-event times for 
all machines regardless of their individual status. To 
determine the next event, the clock module compares the 
smallest of the three failure times with the smallest of 
the three repair times. With this method, there must be a 
means of preventing the program from simulating events con- 


Sidered impossible by the assumptions of the model. 


aah 


For example, one model assumption is that a helicopter 
in cold standby is not subject to failure. It follows that 
the same machine cannot be repaired since it is not in the 
repalr queue. Similarly, a helicopter that is in repair 
cannot fail and a helicopter that is flying cannouwiae 
repaired. To handle these impossible events, the program 
will set the failure time or the repair time, as appropri- 
ate, to equal 1071 which represents 1nimmt,- 

Because of the discrete nature of the graphic display, 
all time values must be converted to integers prior to 
plotting the state-versus-time graph. It is for this rea- 
son that the random number representing the next-event 
time, be it a failure or a repair, is rounded up to the 
next higher integer. The result is that the interarrival 
and service times are not exponentially distributed but are 
instead geometrically distributed. 

The discreteness of the time variable introduces a bias 
in the estimators since they are based on the assumption of 
continuous distributions. The effects of the bias, (ies 
ever, are minimal. The values of the estimates are still 
seen to converge to the limiting values as time advances. 
Despite the generation of geometric random variables, the 
Simulation still provides a sufficient visual demonstration 


of the machine-repairman model. 
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E. RANDOM NUMBER GENERATION 

The author found BASIC's random-number generator, a 
function called RND, unsuitable for simulation. Even so, a 
sample of 10,000 uniform variates produced by RND passed 
with .95 confidence three standard, nonparametric tests for 
uniformity and independence: the serial test, the frequency 
test, and the fone up-and-down test. On the other hand, a 
two-dimensional scatter plot of 5000 pairs of uniform vari- 
ates, shown in Figure 4, clearly illustrate the lattice 


structure of the RND generator. 


U(i-1) 





Figure 4: Two-Dimensional Plot of 5000 Pairs of Uniform 
Variates from the BASIC RND Function. 


MACHREPR uses a random-number generator called RNGEN. 


Written in 8088 assembly language, RNGEN is modification of 
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a generator programmed by Associate Professor B.O. Shubert 
of the Naval Postgraduate School for use on the IBM PC? 
RNGEN uses the algorithm described by Fishman [Ref. 4] 


which is based on the expression 


Mien. = 16807*%. mod UGze ene 
itl a) 


where Xs is the old seed and Xe ay 1s the new seed. 
Dividing the new seed by the modulus produces the uniform 
(0,1), random number. 

This algorithm is the one used by LLRANDOM, a random- 
number generator developed at the Naval Postgraduate School 
and described by Lewis, Goodman, and Miller [Ref. 5] as 
able to produce good quality, uniform variates. 

The simulation computes exponential random numbers by 
the inverse probability integral transformation of uniform 
variates produced by RNGEN. Let U be a uniformly distri- 
buted, random number and let X be exponentially distributed 
with parameter a. An exponential random number is 


generated by setting u = F(x) and solving for x, 


U = PCs ee Bre 
uaa uae 

In (l=u.) =a 

Se iy ee or 


x 


=n (ay oe 
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A. INTRODUCTION 

MACHREPR iS a stochastic, discrete-event simulation of 
a machine-repairman model, which 1S a special case of 
birth-death processes. The model consists of three heli- 
copters, of which two are in service and one is in cold 
standby, with an option of one or two repairmen. The pro- 
gram output is a graphics display composed of a system 
state-versus-time graph, a table of statistics, and animat- 
ed figures that illustrate the dynamics of the process. 

This manual discusses the model and guides the user in 
the step-by-step operation of the program. In describing 
commands and keyboard entries, the following notation is 
used. Single quotes (' ') enclose commands and phrases 
that are to be entered verbatim. The act of pressing a 
Specific key is indicated by brackets (< >) encasing the 
letter or letters on the key. For example, ‘format b:' 
<ENTER> means type the command exactly as it appears 
between the quotes and follow it by pressing the enter key. 
This user's guide shows the commands in the lower case; 


however, the IBM PC will accept upper case letters as well. 


Bee MODEL DESCRIPTION 

MACHREPR contains two versions of the machine-repairman 
model: the one-repairman case and the two-repairman case. 
In both cases, the nominal mission of the model is to keep 


two helicopters flying at all times. If the third 
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helicopter is serviceable, it is kept in cold standby as an 
immediate replacement for the other two. Upon failure of a 
helicopter, a repairman is assigned if he is available; 
otherwise the helicopter joins a queue to await repair. 

Disabled (down) helicopters remain in the repair queue 
until a repairman becomes free. The repair queue ia a 
first in, first out type (FIFO). A newly repaired helicop- 
ter will go immediately into flying operation if less than 
two aircraft are in service. If there are already two 
helicopters flying when a repair is completed, then the 
third goes into the cold-standy queue. 

There are several simplifying assumptions for these 
models. All helicopters are assumed to have exponentially 
distributed service lives with the same failure rate a. If 
T is defined as the length of a helicopter's operating 
life, then the mean time to failure (MTTF) is E(T) = l/a. 
The simulation uses the MTTF as the input parameter that 
establishes the particular exponential distribution for the 
operating life. Helicopters will fly continuously “ane 
failure. A machine that is in cold standby is assumed not 
to be subject to failure. 

Another assumption is that the service time for the 
repairmen haS an exponential distribution with parameter u. 
If S is the time required to complete a repair job, then 
the mean time for repair is E(S) = 1/yu. Like the MTTF, the 


mean time for repair (MTFR) is an input variable for the 
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Simulation. The repairmen are assumed to be equally compe- 
tent and thus have the same MTFR. 

In both cases of the machine~repairman model, the state 
of the system is defined as the number of helicopters down 


Memeime t.. Lhe state variable is N(t) =n, n = 0,1,2,3. 


C. HARDWARE AND SOFTWARE 
1. Hardware 
The simulation is programmed for the IBM Personal 
Computer (IBM PC) with a color monitor. Since the program 
uses color graphics, the microcomputer must be equipped 
with IBM's Color/Graphics Adapter or a suitable graphics 
card from another manufacturer. MACHREPR will not~= run 
without the graphics adapter. The program uses two 
assembly-language subroutines that are loaded in memory 
outside of BASIC's 64 kilobyte (K) workspace; therefore, 
there must be a minimum of 128 K of memory installed in the 
microcomputer. The computer must have at least one 
54-inch, floppy disk-drive. 
2. Software 
MACHREPR and the two subroutines are supplied on a 
distribution diskette. The user must provide the disk 
operating system PC-DOS 2.0 (or 2.10) and Advanced BASIC 
(BASICA) in order to run the simulation. MACHREPR' uses 
BASIC's POKE command and certain program transfers, called 


interrupts, to the Basic Input/Output System (BIOS) in the 
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IBM ROM (read-only memory). Because of these interrupts, 
there can be no assurance that the simulation will operate 
properly on other "IBM compatable" machines. 
3. Program Files 
The distribution diskette contains the following 
program files: 


* MACHREPR.BAS - the main program. 


* RNGEN.SRT - a random-number generator written in 8088 
assembly-language used aS a subroutine to the main 
pEOg Gem. 

* SCRNSHFT.SRT - another assembly-language subroutine 
used to perform a specific task on the graphics 
screen. 


* RNGEN.LST ~- a source-code listing. 
* SCRNSHFT.LST — a souree-code wl istmage 


* MACHREPR.BAT - a batch file that loads and starts the 
Simulation. 


The simulation requires the first three files listed above 
in order to run. The two subroutines are binary-image files 
and can not be executed from PC-DOS. The two source-code 
listings are provided aS a convenience to the uSer. 
MACHREPR.BAT iS a Simple example of how a batch file can be 


used to easily start the simulation. 


D. J*GELTING SlARTED 


1. Making a Backup Copy 


The first step is to make a copy of the distribu- 


tion diskette. After the copy iS made, the original 


38 


diskette should be stored in a safe place and used only to 
make a replacement for an unserviceable application 
diskette. Although there are several ways that the user may 
choose to create an application diskette, this manual will 
Suggest just two: 

* A self-starting (bootable) diskette containing the 
operating system files, the BASICA command file, and 
the simulation's program files. 

* A simple backup copy which contains only the program 
files and which must be used in conjunction with a 
system diskette. 

Owners of microcomputers with a single disk-drive may find 
the first option more convenient. 

To make a bootable diskette, follow the instruc- 
tions in the PC-DOS manual and format a blank diskette 
using the FORMAT command with the /S parameter. The /S 
parameter causes the computer to transfer the operating 
system files from the DOS diskette to the newly formatted 
diskette. The Advanced BASIC file, called BASICA.COM, is 
also found on the DOS diskette and can be copied onto the 
application diskette with the COPY command. Finally, use 
the COPY command again to transfer the simulation's six 
program files from the distribution diskette onto the 
bootable diskette. 

If the second option is desired, then format a 
blank diskette without the /S parameter. This procedure 
does not transfer the operating system files to the appli- 


cation diskette. The next step is to use the DISKCOPY 
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command to make an identical copy of the distribution 
diskette. 
2. Starting the Simulation 

If the first option was used, insert the bootable 
diskette into drive A and turn on the microcomputer. After 
a few seconds, the DOS prompt A> will appear on the screen. 
Tt is here that the batch file MACHREPR.BAT can be used. 
Simply type 'machrepr' <ENTER> to start the simulation. 

If using the non-bootable diskette, then turn on 
the microcomputer with the DOS diskette in drive A. Next 
load Advanced BASIC into memory by typing 'basica' <ENTER>. 
After the computer is in BASIC's command mode, remove the 
DOS diskette from drive A and insert the diskette contain- 
ing the simulation. Type ‘load "machrepr.bas",r' <ENTER> 
to bring the program into memory and to start the 
Simulation. 

If desired, the simulation can be started from a 
second disk-drive. Begin as before by turning on the 
microcomputer with the DOS diskette in drive A and loading 
BASICA. With the application diskette inserted into drive 
B, type ‘load "b:machrepr.bas",r' <ENTER>. During execu- 
tion, MACHREPR reads the program diskette and expects to 
find it in the default disk-drive. When the program fails 
to find the required files on the DOS diskette in drive A, 
it asks the operator for the correct drive label. Simply 


respond to the program's question by pressing <b>. 
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GRAPHIC SIMULATION 
of the 
MACHINE —- REPAIRMAN MODEL 


by R.E. Nelsen 


Submitted in partial fulfillment 
of the requirements for the degree of 


aster of Science in Operations Researc 


from the Naval Postgraduate School 
Monterey, California 


Advisors: J.D. Esary, A.F. Andrus 


Press any key to continue... 





Figure 5: The Title Screen. 


3. Title Screen 
The display will remain dark for a few seconds 
while the program initializes variables, arrays, and func- 
tion definitions. When initialization is done, MACHREPR 


will present the title screen as shown in Figure 5. 


E. THE PROGRAM MENU 

The program menu, depicted in Figure 6, appears next 
and provides the means of controlling the simulation. With 
the options available on the menu, the user may choose to 
read the program instructions, change the model parameters 
or use the default model parameters. The program instruc- 
tions that are available on the screen are a summarized 


version Ore those found in this manual. The model 
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parameters that can be altered by the user are the 
helicopter's mean time to failure, the repairmen's mean 
time for repair, and the number of repairmen to employ. 
Additionally, the Torte lack seed fOr the random-number 


generator can be set. 


PROGRAM MENL 
TYME eer ea Ae Sate See ere aera aan 
lonstructions. 
C>ohange model parameters. 


D>oefault model parameters. 


S>et the random number generator seed. 


E>nd the 


SFeeanieae a ta ee yolk ent Pape eek ae a kok oak ha cadee cdc RL EN 


Enter your selection... 





Figure 6: The Program Menu. 


The simulation starts and ends with the program menu. 
Whenever the user stops the simulation or whenever the sim- 
ulation clock reaches the preset limit of 9950 time units, 
the program returns to the menu. The simulation then can 
be restarted with perhaps a difference set of model parame-. 


ters, or the program can be ended altogether. 
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1. Setting the Random-Number Generator Seed 


To set the seed, press the <s> key. The program 
will then respond by asking for the value of the seed. The 
seed is entered by typing any integer in the range of 1 to 
2,147,483,646 (= 2°+ - 2) and pressing <ENTER>. The pro- 
gram will accept a non-integer value; however, during 
processing the number will be truncated to an integer, with 
the result possibly being different than the desired value. 
After the seed has been entered, the program returns to the 
program menu. 

It 1S not necessary to set the random-number gener- 
ator seed to run the simulation. There is a default seed 
embedded in the generator's program code, and it is updated 
automatically with each call for a random number. If the 
Simulation is stopped and then restarted without setting 
the seed, the generator will use as the initial seed the 
value remaining from the last call in the last simulation 
in . 

2. Changing the Model Parameters 

When <c> 1S pressed, MACHREPR will successively ask 
for the desired values of the mean time to failure, mean 
time for repair, and the number of repairmen. The quanti- 
ties are entered by typing a positive integer and pressing 
<ENTER> in response to each question. For best results, 
the values for the mean time to failure and the mean time 


for repair should be at least ten. The program will accept 
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and function properly with non-integer values of MTTF and 

MTFR; however, these values will be printed on the graphics 

display as integers. Negative numbers will not be accepted 

for the MTTF and the MTFR. The digits 1 and 2 are the vom 

valid responses for the number of repairmen. The program 

will immed tater ys start the simulation after the number of 
repairmen has been entered. 
3. The Default Parameters 

Tf this option is selected, then the simulation 

Starts with the model parameters set at default values. 


These values are: 


MTTF = 15 time units, 
MTFR = 10 time units, 
Number of repairmen = l. 


The default parameters are selected by pressing <d>. 
4. Ending the Program 
To end the program, press <e>. MACHREPR will clear 
the screen, exit the graphics mode, and leave the computer 
in the BASIC command mode. If the user wants to exit BASIC 


and return to the DOS, type 'system' <ENTER>. 


ee THE GRAPHICS pista, 

MACHREPR is designed to give the observer the illusion 
of real-time passage as the probabilistic events occur in 
the machine-repairman system. The graphics display indi- 


cates the model's activities. As shown in Figure 7, the 
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Pane mcs plavetsdivided imtoG three sections: the state- 
versus-time graph, the helicopter display, and the table of 


quantities used with queueing models. 


MACHINE — REPAIRMAN MODEL 
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Figure 7: An Example of the Graphics Display. 


1. The State-Verses-Time Graph 





The upper half of the display contains the state- 
versus-time graph. The values of the state variable N(t), 
the number of machines down at time t, are on the ordinate. 
There are fifty divisions on the abscissa which represent 
units of time. The current time of the system is printed 
memene right-hand end of the abscissa. 

With each unit advancement of the simulation clock, 


the program plots a horizontal line one time unit in length 
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at the height corresponding to the current State or seme 
System. The accumulation of these individual plots forms a 
step graph as shown by the example in Figure 7. When the 
graph reaches the end of the display, the program makes 
room for the next plot by calling SCRNSHFT to™shifttwiaa 
entire graph one unit to the left. In this way, the 
model's activities for the past fifty time units can be 
seen in a single glance. 
2. The Helicopter Display 

The three helicopter images in the lower left 
corner of the graphic display indicate the current state of 
the system. The display in Figure 8 shows how the helicop- 
ters are arranged in the initial state 0. All three 
machines are colored green. Helicopters B and C, seen wima 
the upper half of the display, are flying. Helicopter A on 
the “ground 91s imeodid standby. 

Since all helicopters are initially either flying 
or in cold standby, the next event will be a failure in 
either machine B or C. A machine that fails is grounded 
and transformed into a red colored helicopter with a broken 
rotor. Because there is only one down helicopter, the sys- 
tem is in state 1 and so there is at least one idle 
repairman. The program draws a graphic figure of a man 
beside the down aircraft to represent a machine that is in 
repair. For example, Figure 9 depicts state 1 with 


helicopter B in repair. 
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Figure 8: Initial Helicopter Display in State 0. 
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Figure 9: Helicopter Display in State l. 


As the simulation proceeds, random numbers are gen- 


erated and transformed into interarrival (failure) times 
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and service (repair) times. The program changes’ the 
appropriate graphic figures when the simulation clock 
reaches the time for the next event. 

Figure 10 shows an example of a helicopter display 
that indicates state 2. In this case, C iS in Yepaumae 
Helicopter A is also disabled, and since it does not have a 


repairman, A 1S in the repair queue. 
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Figure 10: Helicopter Display in State 2. 


In Figure ll, all three machines are down, thus the 
system 1S in state 3. There is a repair queue of length 
two in this example. 

It should be noted that Figures 10 and ll are 
illustrations from the one-repairman model since only one 


repairman appears in either display. If the simulation is 
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running the two-repairman model, then the program will draw 


the second repairman beside the helicopter that fails after 


the one already in repair. 
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Figure ll: Helicopter Display in State 3. 


3. The Quantities of Interest 


In the study of queueing systems, there are several 


quantities of interest that are used to evaluate the sys- 


tem's performance. In the lower right corner of 
graphics display (see Figure 7) there is a table of 
more commonly used quantities. 

The row labels have the following definitions: 


* AR - The average arrival rate of machines into 
repair system. 


* Wq - The average waiting time of a machine in 
repair queue. 
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the 


the 


the 


the 


* Lq - The average number of machines in the “repaus 
queue. 


* W - The average amount of time that a machine is in 
the repair system. This quantity is the sum of 
Wq and the expected service time. 

* L - The average number of machines in the repair 
system. This quantity includes those machines 
in the repair queue and those machines in 
‘repair. 

~SAv i= he machine avallability defined as the 
proportion of the time that at least one 
helicopter is flyable. 


* Pd - The probability that a down machine must wait 
for repair. 


* Ps - The probability that an up machine must wait in 
cold standby. 


The numbers that appear in the right-hand column 
under the heading "Limit" are the limiting values of the 
quantities of interest. In other words, the quantities 
converge to the limiting values as time goes to infinity. 
The limiting values are calculated with the balance equa- 


tions that correspond with the particular machine-repairman 


model. 

The values in the left-hand column are estimates of 
the limiting values. These quantities are computed with 
data collected from the simulation. Upon the occurance of 


an event, the program recomputes the estimates after taking 

into account the newly collected data from the last state. 
Because the plots on the state-versus-time graph 

must be given in terms of integer coordinates, all time 


variables must be rounded up to the next higher integer. As 
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a result, the estimators are biased Since they are based on 
the assumption of continuous time when in fact they are 
calculated with data based on discrete increments of time. 
The effects of the bias are minimal; the simulation still 
provides a sufficient demonstration of the convergence of 


the estimates to the limiting values. 


G. KEYBOARD COMMANDS 

There are three keyboard commands available to the user 
as the simulation iS in progress. When <p> (for pause) is 
pressed, the simulation freezes to give the user the oppor- 
tunity to study the graphic display. To continue the simu- 
lation, press <c>. The third command is <s> for stop the 
Simulation. The <s> key will send the program back to the 


program menu. 


H. PROGRAM LIMITATIONS 

MACHREPR always sStartS a simulation under the same con- 
ditions. The initial state is N(0) = 0 with helicopters B 
and C in flight and helicopter A in cold standby. There 
are no provisions in the program for starting the 
Simulation in another state. Restarting the simulation 
from the program menu causes the current time and all data 
collection variables to be set to zero. 

As noted before, all event times are rounded up to the 


next higher integer for plotting the state-versus-time 


Su 


graph. With exponentially distributed interarriva lees 
service times, there iS a non-zero probability that the 
sojourn in any particular state will be less than one time 
unit. The program will still alter the helicopter displ 
properly and compute the estimates correctly; however, the 
State-versus-time graph may display a jump of two states 
because of the program's inability to plot a line of less 
than one time unit in length. The number of times this 
discrepancy occurs can be minimized if the MTTF and the 


MTFR are set at values greater than ten. 
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MACHREPR.BAS FLOWCHART 


Main Program 


START 


| GOSua 39000 | 
E or 
BASICA & 
color/graphics 


Inftialltze 
of ten-used 


var tables 





DEF FN 
EXPONT(U,A) 
DEF FN 
ROUNDUP (A) 


DEF FN 
MIN(A,B,C) 










Display 
Tttie Screen 





Loed random number 
oenerator 


53 











Ofsolay the 
Programs Menu 






Print 
Instructions 
to the screen 


generator seed 
At = i 


Set default for 
ALPHA, 4. NUMRMAN 
Zero varfables 
& arrays 


Compute 
theoret{cal 
statistics 











PRINT 
theoretical 
statfatics 
to screen 
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Oraw graphic 
display and 


ore fn 


arr se 


wy Clock Module 


Set melo A fafl 
time ta ftnfinit 


Compute a*’s 
fafliure t{ne 


Compute C°’s 
faflure time 


Set all repafr 
times to Infinit 


STOSYHELO$ = ‘A°® 


WHILE 
NX TEVNT<9950 













NXTFAIL = 
MIN(FAILTIME(.)) 
NXTREPR = 
MIN(REPRT IME(.)) 


Find whfcnh helo 
willl fafl next: 
FAIL INGHELO$ 


Set NXTFAIL = 
Inftnir 
NXTREPR = 
MINCREPRIIME(.)) 


Find wnfchn helo 
{sa repafred: 
REPAIREOHELOS 


NXTFAIL @ 
MINCFAILTIME(.)) 
Set NXTREPR = 
tnfintt 


Find which peto 
wfll fafl next: 
FAIL INGHELO 


















Find emfcnhn nelo 
{s repafred: 
REPAIREDHELO 











XTFAIL< 
NXTREPR? 


NXTEVYNT © NXTREPR NXTEVNT © NXTFAIL 


Do 


Shift tiae 
grapn to left 
one un{t 


Check keydDoarad 
buf fer 





<P> 


sed ? 
ee 
<S> 
pressed ? 


N 


Ada sojourn time {n state to 
Cumulative tine In state 


LASTEVNTeNXTEYNT 













Upoate 
statistics 
ang print to 
screen 


56 


Clock Module Subrout!nes 


START 
Suor out {ine 
18226 







FAILINGELOF © 
‘ce 


PRINT 
etatistics 
tadle 


Store helo A tn 


Oree STATE = @ 
he lfcooters 





>) 


Faflture Module 





Incr enent 
SUMF ATL 


nmcrement STATE 


= 


Draw & store 


me FIRSTFAILURE «FALSE 
a 
3 


Y 
N 
2 j a 
8 
= 2 


, Faflure Module - Part A 


det] 
U, mu 


AILTIME(AZ) = 1€9312 


Ground fly{ng 2 


REPRT IME (CAS 
+F NE XPONT ( 


~A 


helo & sssfgn 
repafraan 





PSE 


Move standby 
helo to flying 
status 


RE TURN 
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io, Faflure Module - Part 8 


N Y 
Ground flying Ground Melo & 
helo: no ssaign repelrwen 
repairman 


FAILTIME(AS) © 1E31 CALL RNGEN(U) 
REPRTIME(AZ) = 1£3} Ree 


*F NE XPON 


Bet] 
(U,MU 


~m 


FAILTIME(AS) © 1E€31 


N 
neni RETURN 
N 


WAITREPR$ = °8° 





WAITREPRE = °C° 


RE TURN 


* 
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Ney Faflure Module - Part C 


REPRIIME(AS) © 163) 
ATLTIME(AE) © 1£31 







WAITREPRS$ = ‘AS 


WAITREPR$ = 'C° 


RE TURN 
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Fatlure Module Subrout!nes 








PUT(X.105) , GRNHELO 


















Aj e ue ° PUT (X, 143) , REDHELO 
FE XP ONT CU, MU) 
Delay routine: 
count from 1 to 






REPRTIME(2) = Al 
FAILT IME (2)*1€31 







PUT(X, 343) , REDHELO 





PUT(X, 143) , REOHEL OF 


RE TURN 


PuT (xX, 143) , GOWELO 





RE TURN 
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Repafr Module 


Subrout{ne 
3000 
Keg > 
= 'A* 2 





eo Repalr Module - Part A 


ON AZ 







STOBYHELOS © °C 


STOBYMELOS © °A® STOSYHELOS = °8° 


FAILTIME(AZ) © 1€3} 


REPRITME(AZ) © 1633 


| osm 132010 


Che red helo 
to oslo stencby 


STATE © STATE - 3} 


RETURN 
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CALL RNCEN(U) 





{ stare © STATE ~ 1 


RE TURN 


Repafr Module - 


AZ = 3: X = 107 
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Part 8 


Mates 





REPRIIME(AZ) © Time 
@ FMEXPONT(U ba 


STATE @ STATE - ] 


RETURN 


Repafr Module - Part C 


RC 
REPEPRTIME(AIE) = 
TIME * FNEXPCAT(U, MU) 
FATLTIMEC AL) « TIME 
+ FMEXPONT(U, ALPHA) 
Y 
REPRT IME (CAS) = 1631 
N 
ee to os STATE = STATE - 1 
mato ard fly. 
RE TURN 
¥ 
AZ = } AZ e 3 


WAITREPRS$ = °8° 





STATE @ STATE = } 


CALL RNGEN(U) 
RETURN 
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Repafr Module Subrout{nes 






PUT (x, 143) REEL OC? 










LPUT(K 143), GRMMELO | 


RE TURN 





PuT(X, 243) , REDHELC? 


PUT(X,143), GaeHELO 










Delay routine: 
count fro@ i to 
208 


DUT(X,143), GANHELO 


PUT(X.195), GUHELO 


LPUT( 2,143), REDRELO_ 
RETURN 


PUT(X, 143)  REDMELO? 
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APPEND 


MACHREPR.BAS PROGRAM LISTING 


0 ‘MACHREPR.8AS, Yer 3.1, R. E. Nelsen, 30 Aug 1984 

JOCLERRS SAKE W OFF reGES 

2 GOSU8 30000: ‘Check for 8ASICA and color/graphics adaptor 

3 DEFIANT A-Z 

4 DIM P' (3), SUMP! (3) ,GRNHELO (268) ,REDHELO (268) ,REDHELO2 (268) 

5 DIM TIMECHRT (1609) ,MAINSCRN (8002) 

6 DIM FAILTIME! (3) ,REPRTIME! (3) ,STARTQUE (4) 

7 RNGEN=0: TRUE=~1:FALSE=0:U'=0: SCRNSHIFT=384 

8 FMTIS="889. 888° FMT 2S=" 8888": DRVES="A; ° 

9 FIRSTSCRN=TRUEs: WAL TREPRS$="": FAILINGHELOS="":STDBYHELOS="" 

10 REPAIREDHELO$="":S6$ = STRING$(40, 223) 

SQ DEF FNEXPONT!(AL',A2') = - LOB(AL!) /A2! 

U1 DEF FNROUNDUP(A') = INT(A') + 1 

U2 DEF FNMIN'(A',B',C')=(-(A'CB! AND ANCC') EAD F(-(BICA! AND BICC!) aB')+(-(C'CA! AND CHCB!)#C!) 
33 DEF FNHELOPOSIT(A$)=-( INSTR("ABC® ,AS)=1)+(-( INSTR("ABC® A$)=2) #54)+(-( INSTR(*ABC* ,A$)=3) #107) 


60 GOSUB 20000: ‘Title screen 

80 ON ERROR GOTO 30020: ‘Load randoa nuaber generator. 
82 DEF SEG = &HIA00: BLOAD DRVES+*RNGEN. SRT” ,0 

84 ON ERRGR 60TO 0 


90 BLOAD DRVE$+"SCRNSHFT.SRT",384: ‘Load screen shift subroutine 


100 ‘2==2=MAIN PROGRAM MENUsss=s==sssssssssscssssseesesscec= 

120 SCREEN 0,1:COLOR 14,3,0:WIDTH 40:CLS 

121 LOCATE 7,14:PRINT “PROGRAM MENU*:PRINT S6$ 

122 LOCATE 9,1:PRINT "Clnstructions.’ 

125 LOCATE 11,1:PRINT °<C>hange aode! paraseters.” 

124 PRINT:PRINT "<D>efault sodel paraseters." 

125 PRINT:PRINT °<S>et the random nuaber generator seed." 

126 PRINT:PRINT "<E)nd the prograa.*: PRINT S6$ 

130 LOCATE 21,2: PRINT “Enter your selection..." 

ES 

140 GOSUB 10000:AZ= INSTR(*ICDSEicdse",A$) 

141 IF AZ=0 THEN 120 ELSE ON AZ GOTO 20100,200,210,150,900,20100,200,210,150,900 
142 ° 

149 ‘=====RANDOM NUMBER GENERATOR SEED ROUT INE============>== 

150 COLOR 14,0,0:CLS 

151 LOCATE 1,3: PRINT °# SET RANDOM NUMBER GENERATOR SEED #°: PRINT S6$ 

152 LOCATE 4,2:PRINT “Peraissible seed values are integers" 

1535 PRINT ° in the range: 1 to 2147483646." 

154 LOCATE 7,2:1NPUT “Enter the seed value...";A% 

155 IF A# (1 OR A#>2147483446% THEN LOCATE 7,2: PRINT STRING$(39,32): GOTO 154 
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159 DEF SEG = &HLA00 

160 AL'= INT (A8/167772108)sA2'= INT ((AB-AL'€167772108) /69536!) 

142 ‘Poke the seed’s upper 2 bytes into RNGEN’s seed storage. 

163 POKE &H164,Al': POKE &H143,A2! 

164 AB=AB-AL!'#167772108-A2' #45556! sAl'= INT(A8/2546)2A2'= AB-AL'#256 

165 ‘Poke the seed's Tower 2 bytes into RNGEN's seed storage. 

166 POKE &H162,Al': POKE &H141,A2! 

167 6070 120: ‘Return to senu 

168 ° 

199 ‘s=2==PARAMETERS ROUT [NEssssssssssssscssessscssssssszsss 

200 COLOR 15,5,0:CLO:MSG$=")>Enter a positive integer anly.<<* 

202 LOCATE 1,4:PRINT °# CHANGE AGDEL PARAMETERS t°:PRINT S6$ 

203 LOCATE 4,1: INPUT "Enter Mean Tige To Failure...°,A'sIF A'<=0 THEN LOCATE 23,3:PRINT MS6$;:LOCATE 
A, AsPRINT STRINGS (40,32);:6070 203 ELSE ALPHA!=1/A' 

204 LOCATE 6,1: INPUT “Enter Mean Tige For Repair...°,A'sIF A'<=0 THEN LOCATE 23,3:PRINT MS6$;:LOCATE 
6, AsPRINT STRINGS (40,32);:60T0 204 ELSE MU'=1/A! 

205 LOCATE 23,1:PRINT STRINGS(40,32);:MSGS="" 

206 LOCATE 8,1:INPUT “Enter nuaber of repairmen (1 or 2)...°,NUNRMAN: IF (NUMRMAN<>1) AND (NUMRMANK >2 
) THEN LOCATE 8,1:PRINT STRINGS (40,32);:GOTO 206 ELSE 250 

210 ALPHA'=1/15: MU!'=1/10: NUMRMAN=1 

200 NXTFAIL'=O:NXTREPR!=0:F IRSTFAILURE=TRUE: NX TEVNT=0: TIME=0: SUMFAIL=0:LASTEVNT=0: STATE=0: XX=31: SUMW 
B=() 

202 FOR 1=0 TO 3:SUMP! (1)=0:STARTQUE (1)=0:NEXT 

260 ° 


300 P'(1)=2€ALPHA!/MU! 

320 IF NUMRMAN=1 THEN P!(2)=P! (1) €P'(1)sP!(3)=P!(2)#ALPHA!'/MU! ELSE P!(2)=P! (1) ALPHA! /MU!sP! (3) =P! ( 
2) #ALPHA!/ (MU! +AU!) 

ae wOOD=1/ (L4P 8 (1)4P!(2)4P!(3)) 

S2a Pr! ())=P! (1) aP! (0) 

oor ((2)=P'(2)4P! (0) 

327 P!(3)=P! (3) #P! (0) 

330 IF NUMRMAN=1 THEN PDOWN'=P'(2)+P!'(3)3L@'=PDOWN'+P! (3) ELSE L@'=P'(3):PDOWN'=P! (3) 
340 L'=P!(1)+2#P'(2)+34P! (3) 

390 LAMBDA!=(2#P!(0)+2#P!(1)+P!'(2)) €ALPHA! 

360 WE'=LQ!/LAMBDA! 

370 W'=L'/LAMBDA! 

380 PSTDBY'=P! (0) 

390 AVL'=PE(O)H(P (1) 4P! (1) +P! (2))/3 

Sia 

499 ‘ss===PRINT MAIN DISPLAY ROUT NEssssss2esesszeressssez=2 
JOO SCREEN 1,0:COLOR 9,0 

J02 IF FIRSTSCRN THEN 60SUB 15000:FIRSTSCRN=FALSE: 6070 510 
505 PUT (0,0) ,MAINSCRN 

S10 LOCATE 14,27:PRINT USING FHT25;1/ALPHA!; 

511 LOCATE 14,37:PRINT USING FMT2$31/MU!; 

v12 Y=33:605UB 10100; ° Print statistics ta screen. 
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520 LOCATE 12,37:PRINT USING FMT2$3T IME; 

32) GOTO 1000 

eee 

899 ‘s====PROGRAM END ROUT INEssses=ssssesessessccceccccecssss 
900 SCREEN 0,0:WIDTH 80:COLOR 6:PRINT “Program ended... °:END 
ooo 

O09 sasast UC NODC basaaaaeeen ena = ae oa aaa a ee a ee ee 
1000 FAILTIME' (1)=1E +31: “Calculate ist failure tiae (B & C) 
1010 DEF SEG =%HIAOO:CALL RNGEN(U!) 

1012 FAILTIME'(2)= FNEXPONT! (U! ALPHA!) +TIME 

1020 CALL RNGEN(U!) 

1022 FAILTIME! (3)= FNEXPONT! (U! ,ALPHA!) +TIME 

1030 REPRIIME! (1)=1E+31: “LE+31 = “infinity” 

1032 REPRTIME! (2)=1E+31 

1034 REPRTIME'(3)=1E+31 

1034 STDBYHELOS="A° 

1099 ‘## Deteraine next event ## 

1100 WHILE NXTEVNT<9950; ‘Sets upper limit to prow run tiae. 
1110 IF NOT STATE=0 THEN 1120 

1112 NXTFAIL!= FNMIN'(FAILTIME' (1) FALLTIME! (2) ,FAILTIME! (3)) 
1114 NXTREPR! = 1E+3! 

1116 GOSUB 10210:G0TO 1180 

1120 IF NOT STATE=3 THEN 1130 

1122 NXTFAIL'= 1E+31 

1124 NXTREPR'= FNMIN' (REPRTIME!' (1) ,REPRTIME!' (2) ,REPRTIME! (3)) 
1126 GOSUB 10220:60T0 1180 

1130 NXTFAIL'= FNMIN'(FAILTIME! (1) PAILTIME!' (2) FAILTIME! (3)) 
1132 NXTREPR'= FNMIN' (REPRTIME!' (1) ,REPRTIME! (2) ,REPRTIME! (3)) 
1134 GOSUB 10210 

1136 GOSUB 10220 

1180 [F NXTFAIL'<=NXTREPR! THEN 1184 

1182 NXTEVNT = FNROUNDUP(NXTREPR'):60TO 1200 

1184 NXTEVNT = FNROUNDUP(NXTFAIL!) 

plese. 

1199 *## Run clock and update time chart ## 

1200 DEF SEG:POKE &HbA,0: “Clear key buffer 

1210 WHILE TIME ¢ NXTEVNT 

1220 TIME=TIME+! 

1230 1F TIME<=50 THEN GOSUB 11050: XX=XX+5:FOR [=1 TO 9OO:NEXT:G60TO 1240 
1232 DEF SEG = &H1A00 

1234 CALL SCRNSHIFT 

1238 XX=2746:GOSUB 11050 

1240 LOCATE 12,37:PRINT USING FMT2$; TIME: 

1241 AS=INKEYS sIF AS="* THEN 1240 

1242 1F A$="s*° OR A$="S* THEN 120 

1244 1F A$="p* OR AS="P* THEN GOSUB 10010 

1250 DEF SEG:POKE &H6A,0: ‘Clear key buffer 

1260 WEND 

P76 ly 

1299 ‘## Update display and determine next event ## 

1300 SOUND 120,3:SUNP! (STATE) =SUMP! (STATE) +1 1ME-LASTEVNT 
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1302 LASTEVNT = NXTEVNT 

1310 IF NXTFAIL'<=NXTREPR! THEN GOSUB 2000 ELSE 60SUB 3000 

foil 

1399 ‘t# Update statistical estiaates t+ 

1400 LAMBDA'=SUMFAIL/TIME 

1410 IF NUMRMAN=1 THEN PDOWN'=(SUMP! (2)+SUMP!(3))/TIME:LQ!'=PDOWN'+SUMP! (3)/TIME ELSE LQ'=SUMP! (3)/TI 
ME: PDOWN'=L@! 

1420 L'={SUMP! (1) +SUMP! (2)4SUMP! (2) +SUMP! (3) +SUNP! (3)+SUMP!(3))/TIME 

1430 WO'=SUMWQ/SUNFAIL: W!=WO!' +1 /MU!': PSTDBY!=SUMP! (0) /TIME 

1440 AVL'=SUMP! (0) /TIME+(SUNP! (1) +SUMP! (1) +SUMP!(2))/(TIME+TIME+TIME) 

1450 Y=25:60SUB 10100: ‘Print estig@ates to screen 

1451 ° 

1490 A$= INKEYS$:1F A$="° THEN 1500 

1491 IF Ag=*s* OR A$="S" THEN 120 

1493 IF Af="p" OR A$="P" THEN GOSUB 10010 

1495 DEF SEG: POKE &HSA,0: ‘Clear key buffer 

1500 WEND :‘No stop command so go back to line 1100. 

1510 CLS: LOCATE 12,1: PRINT “Program time li@it reached.* 

1512 PRINT "Press any key to return to menu...* 

1514 GOSUB 16000: GOTO 120 

5) 

ame nen leeas=s—e—=a—e—=s——esnesessoseesste==s= 

2000 SUNFAIL=SUMFAIL+L:STATE = STATE + ! 

2020 IF FIRSTFAILURE THEN GOSUB {2100:FIRSTFAILURE=FALSE:60T0 2130 

2030 AL= INSTR{*ABC* ,FAILINGHELOS) 
2032 X= FNHELOPOSIT(FAILINGHELOS) 
2040 ON STATE GOTO 2100,2200,2300 
2041 ° 

2100 DEF SEG =8HLAOO:CALL RNGEN(U!) 
2110 REPRTIME! (AZ)= FNEXPONT! (U! MU!) +T IME 
2412 FAILTIME!' (AZ)= 1E+31 

2120 GOSUB 12220: ‘Ground flying helo & assign repairaan. 
2130 Ad= INSTR(*ABC* ,STDBYHELOS) 

2132 X= FNHELOPOSIT (STDBYHELOS) 

2140 CALL RNGEN(U!) 

2142 FAILTIME! (AZ)= FNEXPONT! (U! ,ALPHA!)+TIME 

2150 GOSUB 12230: ° Move stand-by helo to flying status. 
2140 STOBYHELO$="°s RETURN 

Zils” 

2200 IF NUMRMAN =2 THEN 2240 

2219 60SUB 12210: “Ground flying helo: no repairman. 
2220 WAITREPRS = MIDS(*ABC",AZ,1) 

2222 REPRTIME! (AZ) = LE+3isFAILTIME! (AZ) = 1E+31 
2224 STARTQUE(AZ) = TIME 

2230 RETURN 

2240 6O0SUB 12220: ‘Ground flying helo & assign a repairean. 
2250 DEF SEG =%HIA00: CALL RNGEN(U!) 

2252 REPRTIME! (AZ) = FNEXPONT! (U!,MU!)+TIME 

2254 FAILTIME! (AZ) = LE+31: WAITREPRS="" 

2260 RETURN 

ezol 


wae 


2500 
2310 
Polo 
2320 
2330 
25a 


BIo8e 


5000 
5010 
5020 
5021 
5100 
3110 
3120 
5130 
5140 
S141 
3200 
3210 
3212 
3220 
3230 
3240 
5242 
3244 
$250 
5260 
5262 
3270 


REPRTIME' (AZ) = JE+31:FAILTIME' (AZ) = 1£+31 

IF NUMRMAN=2 THEN WAITREPRS= MIDS(“ABC*,AZ,1) 
STARTQUE(AX) = TIME 

GOSUB 12210: ‘Ground flying helo: no repair@an. 
RETURN 


A= INSTR(*ABC* ,REPAIREDHELOS) 
X = FRHELOPOSIT (REPAIREDHELOS) 
ON STATE GOTO 3100,3200,3300 


STDBYHELOS= MID$(°ABC*,AZ,1) 

FAILTIME! (AZ) = 1E+31 

REPRIIME'{AZ) = 1&+31 

G0SUB 13010: ‘Change red helo to cold stand-by helo. 
STATE = STATE-i: RETURN 


DEF SEB =8HIAOOsCALL RNGEN(U!) 

FAILTIME! (AZ) = FNEXPONT' (U!, ALPHA!) +TIME 
REPRIIME'(AZ) = LE+t31 

GOSUB 13020: ‘Change red helo to green and fly. 

IF NUMRMAN=2 THEN STATE=STATE-1: RETURN 

AZ= INSTR(*ABC* WAI TREPRS) 

X = FNHELOPOSIT(WAITREPRS) 

SUMWQ = SUMWQ+TIME-STARTQUE(AZ) 

GOSUB 13030: ‘Change red helo to red helo w/ repairaan. 
CALL RNGEN(U')sREPRTIME! CAZ)= FNEXPONT! (U! MU!) +T IME 
WAI TREPRS=°° 

STATE=STATE-1: RETURN 


ole 


3300 
5310 
$312 
3520 
$330 


Iti 
JJd 


3334 
3340 
3342 
3344 
3346 
3330 
3360 
3362 
3364 
3370 
3372 
3374 
3380 
3382 
3390 


DEF SEG =&HIAOO:CALL RNGEN(U!) 

FAILTIME!' (AZ)= FNEXPONT! (U! ALPHA!) +TIME 

REPRTIME! (AZ)= 1E+31 

GOSUB 13020: ‘Change to green helo and fly. 

Aix= INSTR(*ABC® ,WAITREPRS) 

X= FNHELOPOSIT(WAITREPR$): ‘Helo at head of queue. 
SUMWO = SUMWO+TIME-STARTQUE(AIZ) 

GOSUB 13030: ‘Change red helo to red helo w/ repairaan. 
CALL RNGEN(U!) 

REPRTIME! (ALZ)= FNEXPONT'(U! MU!) +TIME 

IF NUMRMAN=2 THEN WAITREPR$="*°:60T0 3390 

ON AR GOTO 3362,5372,3382 :° Advance the queue. 
"REPAIREDHELO$ = °A° 

IF WAITREPRS="B° THEN WAITREPRS="C° ELSE WAITREPRS$="B" 
6070 3390 

"REPAIREDHELO$ = °B* 

IF WAITREPRS="A* THEN WAITREPRS="C* ELSE WAITREPRS$=°A° 
GOTO 3390 

"REPAIREDHELOS = °C* 

IF WAITREPR$="A* THEN WAITREPRS$="B° ELSE WAITREPRS=°A° 
STATE=STATE-1: RETURN 
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Sook 

9998 ‘stst# SUBROUTINES FEEEEEEESEESEESEEESERESELELEL EL EEE SE 
a999 

10000 DEF SEG:POKE &H4A,0: ‘General input routines=s=ssss-=s= 
10001 Ag= INKEY$:IF AS="" THEN 10001 ELSE RETURN 

10002 ° 

10010 DEF SEG: POKE &H4A,0: ‘Pause routines==sssssssssssssss= 
10011 AS= INKEY$: IF A$="c" OR A$="C" THEN RETURN ELSE 60TO 10011 
1012 © 

10099 ‘=====STATISTICS PRINT SUBROUTINE========s===ss=ss=s2== 
10100 LOCATE 18,Y:PRINT USING FAT1$;LANBDA!; 

10101 LOCATE 19,Y:PRINT USING FATI$;WQ!; 

10102 LOCATE 20,Y¥:PRINT USING FMT1$;L@!: 

10103 LOCATE 21,Y:FRINT USING FATI$;:W!; 

10104 LOCATE 22,Y¥:PRINT USING FMTI$;L!; 

10105 LOCATE 23,¥:PRINT USING FATI$;AVL!; 

10106 LOCATE 24,Y¥:PRINT USING FAT1$;PDOWN'; 

10107 LOCATE 25,Y:PRINT USING FATIS;PSTDBY!; 

10108 RETURN 

10109 ° 


10210 IF NXTFAIL'=FAILTIME! (1) THEN FAILINGHELOS="A° ELSE IF NXTFAIL‘=FAILTIME! (2) THEN FAILINGHELOS 
="B° ELSE FAILINGHELOS="C° 

10211 RETURN 

10220 IF NXTREPR'=REPRTINE! (1) THEN REPAIREDHELOS="A" ELSE IF NXTREPR'=REPRTIME!'(2) THEN REPAIREDHEL 
O$="B° ELSE REPAIREDHELOS="C* 

10221 RETURN 

10222 ° 

ee neem None SIbK Nb ss]s==-2=sososs=S==ncSss=s>s> 

11050 AZ=STATE+1:ON AZ GOTO 11051,11052,11053,11054 

11051 Y=75:LINE (Xx,Y) - STEP (5,0),2:RETURK 

Mee Y=9:LINE (XX,¥) — STEP (5,0) ,2: RETURN 

11053 Y=43:LINE (X1,¥) - STEP (5,0),2:RETURN 

11054 Y=27:LINE (XX,Y) - STEP (5,0) ,2:RETURN 

11055 ° 

12099 ‘s===s=FIRST FAILURE SUBROUTI  NEs2s2s2sss2sse-ssccserze2 

12100 DEF SEG =&HIAQO:CALL RNGEN(U!):A'= FNEXPONT'(U! MU!) +TIME 

12102 IF FAILINGHELO$="B" THEN X=54:X1=7B:REPRTIME! (2)=A!'sFAILTIME! (2)=1E4+31:60T0 12106 
12104 IF FAILINGHELOS="C° THEN X=107:X1=151:REPRTIME! (3)=A'sFAILTIMNE! (3)=1E+31 

12106 FUT (X,105) ,GRNHELO 

12110 PI = 3.141593 

Weve CIRCLE (i! ,161),12,2,0,P1 

12112 DRAW °C2 BUI2 U4 M-12,+10 M-8,-4 DIO BM#20,-16 M#12,+10 Rb M+4,+5° 

12113 CIRCLE (11,165) ,14,2,7.5#P1/6, 11. 54P1/6 

12114 DRAW "C2 BD2 L4 NM-2,-15 M-8,-A":DRAW "D11 M-4,4+5 D4 NLJ NRI U3 M#4,-3° 

12115 DRAW °BR20 M+4,+3 D3 NLI MRI U4 M-4,-5 ULI": DRAW °M-B,+6 NM+2,-15 L4" 
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LETTE PRINT STEP AO, alae 

[2117 DRAW °CO L4 DS Ri Ub Rb DS RI US L4” 

12113 CIRCLE STEP 4G 42,0; oblecler WCU) ae 

12120 GET (X,143)- STEP (48,40) ,REDHELO: X1=X1+19 

12130 CIRCLE (x1, 160) 4 PPR INC Togo 

W213) DRAW. °C] BA=$,-2 Nugent ns 2 

12132 LINE. STEP (0;6}= STEP (2al ina. be 

12133 LINE STEP (1,0)- STEP (-4,6),1,BF 

12134 LINE STEP (120)-" SUER mt? 6) ce 

12135 DRAW “CI L3 BUY LZ R-4,-1 UL Ata tl R2°SPAINT SEEP elie 
12140 GET (X,143)- STEP (4B, 40) ,REDHELO2:RETURN 

12141 ° 

12209 =====FAILURE DRAWING SUBROUTINES===2s==2=ss=2=2s2s2s2s== 

12210 PUT (X,105) ,GRNHELO:PUT (X,145) ,REDHELO:RETURN: ‘ No repairean 
Petia 

12220 PUT (X,105) ,GRNHELO:PUT (X, 145) ,REDHELO 

12221 FOR 1=1 TO SOO:NEXT:PUT (X,143) ,REDHELO 

12222 PUT (X,143) ,REDHELO2:RETURN: © With repairaan 

122234 

12230 PUT (X,143) ,GRNHELO:PUT (X,105) ,GRNHELO:RETURN: ° Cold standby to flying status. 
12231 ° 


13009 ‘=====REPAIR MODULE SUBROUTINES==:=ss:s2s2sssssssesee= 
13010 PUT (x, 143) REDHELO2:PUT (X,143) ,GRNHELO:RETURN: ° To cold standby 
13011 


13020 PUT (X,143) ,REDHELO2Z: FUT (X,143),GRNHELO 

13021 FOR 1=1 TO 200:NEXT: PUT (X,143) ,GRNHELO 

13022 PUT (X,105),GRNHELO: RETURN: ‘Repaired and goes flying. 

12023 

13030 FOR J=1 TO 200:NEXTsPUT (X,{43) ,REDHELO: PUT (X,143),REDHELO2: RETURN: “Assign a repairaan. 
13031 ° 

14999 ‘s====MAIN SCREEN DRAWING SUBROUTINE=======ss=2s2s==== 

15000 LINE (0,8)-(319,8),2:LINE (0,10)-(319,10),2 

{5002 LINE (31,23)-(31,80):FOR Y=27 10 73 STEP Lo:LINEIZ7 ysis eee 

15004 LINE (31,80)-(281,890) 

15006 FOR x=36 TO 28! STEP § 

15008 IF (X=81) OR (X=131) OR (X=1B1) OR (X=231) OR (X=281) THEN LINE (X,80)-(X,85) ELSE LINE (x,80) 
-(X,83) 

15010 NEXT 

15012 LINE (0,100)-(319, 100) ,2:LINE (0, 102)—(Sio loz 

15014 LINE (161,102)-(161,199),2 

15016 LOCATE 1,8:PRINT *MACHINE - REPAIRMAN MODEL’ 

15018 LOCATE 4,3:PRINT °S°sLOCATE &6,3:PRINT °2°:LOCATE 8,3:PRINT “1°sLOCATE 10,3:PRINT °0° 
15020 LOCATE S,i:PRINT °N(t)*sLOCATE 12,32:PRINT “Tise:°; 

15022 ‘Helocopter drawing routine 

15024 Pl = 3.141593: %=25: Y2=165 

15026 6OSUB 15038 

15028 X=78:¥2=127 
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15030 GOSUB 15038 

15032 1=131 

15034 GOSUB 15038 

15036 GOTO 15062 

15038 Y1=Y2-4 

15040 CIRCLE (1,Y1),12,1,0,PI 

15042 DRAW "C1 BUI2 U4 L22 R44" 

15044 CIRCLE (X,¥2),14,1,7.5#P1/6, 11, 5#P1/4 

15046 DRAW "C1 BD2 L4 NN-2,-15 H-8,-4° 

15048 DRAW "D11 M-4,45 D4 NLJ NRI U3 M+4,-3° 

15050 DRAW *BR20 M+4,+3 D3 NL1 NRI U4 M-4,-5 ULL" 
15052 DRAW *M-8,+6 NMN#2,-15 L4" 

15054 PAINT STEP (0,4) ,1,! 

15056 DRAW "CO L4 Do Ri Ud RA DO RI Ud LA 

$5058 CIRCLE STEP (0,4),1,0:PSET STEP (0,0) ,0 

15060 RETURN 

15062 ‘Statistics Routine 

15064 LOCATE 14,22:PRINT "MITF:";:LOCATE 14,32:PRINT “MTFR:*; 
15066 LINE (161,113)-(319, 113) ,2 

15068 LOCATE 16,22:PRINT *Var*:LOCATE 14,24:PRINT “Estiaate’:LOCATE 16,35:PRINT “Liait" 
15070 LOCATE 18,22:PRINT "AR" 

15072 LOCATE 19,22:PRINT "Wg" 

15074 LOCATE 20,22:PRINT "Lq° 

15076 LOCATE 21,22:PRINT "W" 

15078 LOCATE 22,22:PRINT "L" 

15080 LOCATE 23,22:PRINT "Avi" 

15082 LOCATE 24,4:PRINT "A": :LOCATE 24,10:PRINT *B": 
15084 LOCATE 24,17:PRINT *C*s:LOCATE 24,22:PRINT "Pd°: 
15086 LOCATE 25,22:PRINT "Ps"; 

15087 GET (0,0) - (319,199) ,MAINSCRN 

15088 GET (1,143)-(49, 183) ,6RNHELO 

15089 RETURN 

15090 


20000 SCREEN 0,1: COLOR 14,1,1: WIDTH 40: CLS 

20010 LOCATE 2,12: PRINT "GRAPHIC SIMULATION" 

20012 LOCATE 4,18: PRINT “of the" 

20014 LOCATE 6,9 : PRINT “MACHINE - REPAIRMAN MODEL" 

20016 LOCATE 9,14: PRINT “by R.E. Nelsen" 

20020 LOCATE 14,5: PRINT "Subaitted in partial fulfillaent* 

20021 LOCATE 15,2: PRINT “of the requiresents for the degree of" 
20022 LOCATE 17,1: PRINT "Master of Science in Operations Research’ 
20023 LOCATE 19,4: PRINT "from the Naval Postgraduate School" 

20024 LOCATE 20,11: PRINT "Monterey, California’ 

20025 LOCATE 22,2: PRINT “Advisors: J.D. Esary, A.F. Andrus” 

20030 LOCATE 25,2: PRINT “Press any key to continue...";:60SUB 10000 
20040 CLS:LOCATE 1,14: PRINT “Acknowledgaent* 

20041 LOCATE 3,4:PRINT "The author would like to thank" 

20042 LOCATE 7,2:PRINT “Associate Professor Bruno Shubert for” 
20043 LOCATE 9,2:PRINT "the donation of his asseably-Ianguage" 


u 


20044 LOCATE 11,2:PRINT ‘randoa-nuaber generator which foraed’ 

20045 LOCATE 13,2:PRINT “the basis for RHGEN.SRT.* 

20044 LOCATE 25,2:PRINT "Press any key to continue...°;:G0SUB 10000: RETURN 

20047 ° 

20099 ‘ss=s=INSTRUCT I ONS=s22=2ss22sssssssssscsssscsccsssssssz 

20100 SCREEN 0,1:COLOR 14,1, 1:WIDTH 40;CL5 

20110 LOCATE 1,4:PRINT “Instructions for the Prograe Menu’; 

20111 LOCATE 3,2:PRINT “Press 4s) to set the randoe-nuaber seed*;:PRINT “and enter a positive intege 
r when asked. The prograa returns to the aenu." 

20112 PRINT:PRINT ° Press (c> to change aodel parameters. The prograa will ask for the aachine’s a 
ean time to failure (MTTF), the repair-°; 

20113 PRINT “men's aean tiae for repair (MTFR), and the nusber of repairaen to eaploy. The" 

20114 PRINT °MTTF and the MTFR should be positive integers that are greater than 10. The nuaber 

of repairaen can be | or 2. The simulation starts after the nuaber of repairaen has been entered. 


20115 LOCATE 18,2: PRINT "Press (d> to use default parameters. ";:PRINT "The siaulation will itemediat 
ely start with the values: MTTF = 15, ATFR = 10, nuaber of repairaen = 1." 

20116 LOCATE 23,2: PRINT *To end the prograa, press (e>.° 

20117 LOCATE 25,2: PRINT "Press any key to continue...";:G0SUB 10000 

20120 CLS:LOCATE 1,13: PRINT “Keyboard Commands’; 

20121 LOCATE 4,2: PRINT “During the sigulation, you aay use the’s;:PRINT “following keyboard coamands 


20122 LOCATE 7,5: PRINT “<p> - pause the simulation.’ 

20123 LOCATE 9,5: PRINT "<c> - continue the siaulation.” 

20124 LOCATE 11,5: PRINT °<s> - stop the siaulation and return’ 
20125 LOCATE 13,11: PRINT *to the program aenu.° 

20126 LOCATE 25,2:PRINT “Press any key to continue...";:G0SUB 10000 
20130 CLS:LOCATE 1,13: PRINT “Graphic Display’; 

20131 LOCATE 3,2:PRINT “The variable names in the display are*s:PRINT “defined as:°; 
20132 LOCATE 6,1: PRINT "TIME: current tise on siaulation clock.’; 
20133 LOCATE 8,1: PRINT "Nit): no. of machines down at tise t.°; 
20134 LOCATE 10,3: PRINT "AR: average arrival rate of machines’; 
20135 LOCATE 11,7: PRINT "into the repair queue."; 

20136 LOCATE 13,3: PRINT "Wg: average waiting time of machines’; 
20137 LOCATE 14,7: PRINT “in the repair queue.’s 

20138 LOCATE 16,3: PRINT “La: average length of repair queue.*; 
20139 LOCATE 18,4: PRINT "W: average down tiae of a machine."; 
20140 LOCATE 20,4: PRINT “L: average no. of down aachines.’; 

20141 LOCATE 22,2: PRINT “Avil: machine availability.°: 

20142 LOCATE 25,2: PRINT "Press any key to continue...°;:G0SUR 10000 
20143 CLS:LOCATE 1,3: PRINT "Pd: probability that a down machine’; 
20144 LOCATE 2,7: PRINT “aust wait for repair.°; 

20145 LOCATE 4,3: PRINT °Ps: probilility that an up aachine’; 

20146 LOCATE 5,7: PRINT “aust wait in cold standby.’; 
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20147 LOCATE 8,1: PRINT ° The nugbers listed under the coluan heading ‘Liait’ are the ligiting va 
lues as tiae goes to infinity of the above defined quantities. These values are computed with t 
he balance equations for the saodel.* 

20148 LOCATE 15,1:PRINT ° The values listed under the coluan heading ‘Estimate’ are the estirate 
s of the limiting values. These quantities are calculated with data collected froa’ 

20149 PRINT “the sigulation. The estimates are recoa-puted upon every occurance of an event." 

20150 LOCATE 25,2:PRINT “Press any key to return to menu... °3:GOSUB 10000:G0T0 120 

20151 ° 

29999 ‘=====BASICA AND COLOR/GRAPHICS ADAPTER CHECK========= 

30000 DEF SEG =O0:1F {PEEK(%H410) AND &H3Q) <> &HSO THEN DEF SEG:60T0 30007 

30003 LOCATE 3,1:PRINT “Sorry...” 

30004 PRINT "You do not have the color/graphics sonitor adapter'* 

30005 PRINT "This siaulation uses graphics and requires that adapter." 

30005 DEF SEG : END 

30007 ON ERROR GOTO 30008: PLAY "pid*:6070 30011 

30008 WIDTH 80:CLS:LOCATE 3,1 

30009 PRINT “This siaulation uses advanced BASIC.* 

30010 PRINT “Reload this prograa after using the command ‘BASICA’.°:END 

30011 ON ERROR GOTO 0: RETURN 

30012 ° 


30020 IF ERR <> 53 THEN GN ERROR GOTO 0 
30022 CLS:LOCATE 3,1:PRINT "Qn which drive can this programs be found?":60SUB 10000 
30024 DRVES=AS+"s": RESUME 82 
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© RNGEN.SRT -- BASIC USR Routine 


by B.Q. Shubert 


; Modified by R.E. Nelsen 
; 2 May 1984 


* This subroutine generates and returns a unifore 

; rando@-nuaber in BASIC’s single-precision format. 
; The seed X 1s not iaported nor returned to BASIC, 
; Dut 1t 15 autoraatically updated. The algoriths 

; used 15: 


XY = AY nod a 


; with A = 7*5 and M = 2°31 - 1. RNGEN is called 


froa BASIC using the statesents: 
RNGEN = 0 
DEF SEG = &HiA00 
CALL RNGEN (U!) 


© U' is the returned randp# nuaber and sust be 
; defined prior to the first call to RNGEN. 


SMUSTITUI TITAS TTR TTT TUT UTS T TTT asa a aaa 


3 43 


CSEG = SEGMENT 
ASSUME CS:CSEG,DS:CSEG,ES:NOTHING 


ORG = 0) 
RNGEN PROC =FAR 

PUSH BP ssave for BASIC 
a HOV = BP, SP spoint to args on stack 

PUSH DS ssave for BASIC 
C8 MOV bi-CS smake data accessable 
D8 MOV =s-«&S, AX 

J 

O161 R MOV = AX,LOSEED = =©sget old X 
16 0163 R MOV «= DX, HISEED 
FQ MOV =o ST, AX 
FA MOV =—soDI, DX 
FA 01 CMP sé, 1 scheck if ASX ¢ M 
bA JA MULT 48 348 bit multiply if not 
05 JB MULT32 selse 32 bit aultiply 
F31D CMP =s AX, OF S1DH 

JNB  MULT48 
EO MULTS2; SAL Ak smultiply X by 7*5 
D2 REL eee gt 
EQ SALA, I 
D2 RCL =—sDX, 1 
EQ SAL Aka! 
D2 RC! DX, 
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0025 
0020 
002F 
0031 
0033 
0035 
0037 
0039 
0038 
003D 
O03F 
0041 
0043 
0045 
0047 
0049 
0048 
004D 
OO4F 
0051 
0053 
0055 
0037 
00359 
0058 
0050 
O05F 
0041 
0063 
0049 
0067 
0089 
0048 
004D 
O04F 
0071 
0073 
0075 
0077 
0079 
0078 
007D 
007F 


0082 
0084 
0084 
0088 
0084 


2B 
1B 
8B 
8B 
DI 
Di 
DI 
Di 
Di 
DI 
28 
1B 
GB 
OB 
Di 
Di 
DI 
DI 
Dt 
Di 
28 
1B 
GB 
8B 
DI 
D1 
Di 
Di 
D1 
DI 
28 
1B 
GB 
OB 
Di 
DI 
DY 
D1 
Dt 
Di 
28 
1B 


EQ OL2F R 


Cé 
D7 
FQ 
FA 
E) 
D2 
E) 
D2 
EO 
D2 
C4 
D7 
FO 
FA 
EO 


a 
£ 


EO 


n 
é 


EQ 
D2 
Cé 
D? 
FQ 
FA 
E) 
D2 
E) 
D2 
EO 
D2 
Cé 
D7 
FQ 
FR 
E0 
D2 
E0 
D2 
EO 
D2 
Cé 
D? 


33 DB 


8B 
Di 
Di 
Di 


CB 
EO 
D2 
D3 


MULT48: 


SUE 
SBE 
MOV 
OV 
SAL 
RCL 
SAL 
RCL 
SAL 


RCL - 


SUB 
SBB 
MOV 
OV 
SAL 
RCL 
SAL 
RCL 
SAL 
RCL 
SUB 
SBR 
OV 
OV 
SAL 
RCL 
SAL 
RCL 
SAL 
RCL 
SUB 
SBB 
HOV 
OV 
SAL 
RCL 
SAL 
RCL 
SAL 
RCL 
SUB 
SBB 
JEP 


XOR 
HOV 
SAL 
RCL 
RCL 


AX, SI] 
DX, DI 
S1, Ax 
DI, DX 
AX,1 
DX, 1 
AX,1 
DX, 1 
AX,1 
DX, 1 
AX,S] 
DX, DI 
51, Ax 
DI, DX 
AX, 1 
DX,1 
AX, 1 
DX, 1 
AX, 1 
DX,1 
AX, SI 
DX, DI 
S1, AX 
DI, DX 
AX, 1 
DX, 1 
AX,1 
DX, 1 
AX, 
DX, 1 
AX, 51 
DX, DI 
SI, Ax 
DI, DX 
AX, 1 
DX, 1 
AX,1 
DX! 
AX, 1 
DX, 
AX, SI 
DX, DI 
UNIFO! sdone since AgX<M 


BX, BX seultiply X by 7*5 
CX, BY 

AX,1 sby multiplying X by 8 
DX, 

BX, 


Vo 


908C 
O08E 
0090 
0092 
0094 
0096 
0098 
0094 
009C 
009E 
00A0 
N0AZ 
00A4 
D0A6 
00A8 
NOAA 
OOAC 
NOAE 
0080 
ODOR? 
00B4 
0086 
0088 
NOBA 
OORC 
OOBE 
NOLO 
0002 
00C4 
0004 
N0C§ 
NOTA 
00cC 
OOCE 
90D0 
99D2 
00D4 
00D 
0008 
OODA 
DC 
QNDE 
O0E0 
ONE2 
OOE4 
00E6 
NOER 
OOEA 
NOEL 


DI 
DI 
DI 
Dt 
Di 
Di 
2B 
1B 
1B 
8B 
8B 
8B 
DI 
DI 
D1 
D1 
D1 
Di 
D1 
Di 
Di 
2B 
1B 
18 
9B 
8B 
8B 
DI 
Dl 
Di 
D1 
D1 
DI 
DI 
Dt 
Di 
28 
1B 
1B 
8B 


E) 
D2 
DS 
EQ 
D2 
D3 
Cé 
D/ 
D9 
Ff) 
FA 
C8 
E0 
D2 
DS 
EQ 
D2 
D3 
E0 
D2 
D3 
Cb 
D? 
D9 
FQ 
FA 
CB 
ED 
D2 
D3 
EO 
D2 


5 
3 


E) 
D2 
D3 
Cé 
D/ 
D9 
FQ 
FA 


BE CR 


Di 
Dt 
DI 
Dt 
D1 
DI 
DI 


EO 
D2 
D3 
EQ 


? 
DS 
E0 


SAL 
RCL 
RCL 
Sé 

RCL 
REL 
SUB 
SBB 
OBB 
MOV 
MOV 
OV 
SAL 
RCL 
REL 
SAL 
RCL 
RCL 
SAL 
REL 
RCL 
SUB 
SBB 
SBB 
NOV 
MOV 
MOV 
SAL 
REL 
RCL 
SAL 
RCL 
RCL 
SAL 
RCL 
RCL 
SUB 
OBB 
SBB 
MOV 
MOV 
MOV 
SAL 
RCL 
RCL 
SAL 
REL 


SAL 


AX, ! 
Dx, 
BX, ! 
aX, 1 
DX, 1 
BX, 1 
aX, SI 
DX,DI 
BX, CX 
S1, AX 
D1, DX 
CX, BY 
AY, 1 
DX, 1 
BY, 1! 
AX, 1 
DX, 
BX, ! 
aX, ! 
DX, 1 
BX, 1! 
AX,S1 
DX,DI 
BY, CY 
SI, AY 
DI, DY 
CX, BY 
aX, ! 
Dx, ! 
BX, ! 
AY, 1 
Di, 1 
BX, 1! 
AX, ! 
DX, 
BX, 1 
AX,SI 
DY, DI 
BX, CX 
S1, AX 
DI, DY 
CX, BY 
aX! 
DX, ! 
BX, 1! 
AX, 1! 
DX! 
Bi, 
Ax, 1 
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;subtract original X 


iSaving new X 


vand repeat this & tises 


OOEE 
0OF 0 
O0F2 
OOF 4 
NOF 6 
OOFB 
OOFA 
OOFT 
OOFE 
0106 
0102 
0104 
0104 
0108 
010A 
O10C 
O19E 
0110 
0112 
0114 
0118 
0118 
OLLA 
O11C 
OLLE 
0121 
eS 
0127 


012A 
0128 


012F 
0132 
0134 
0138 
013A 
013C 
0136 
0140 
0142 
0144 
0185 
0148 
0148 
O14C 


D1 D2 

D1 D3 

2B Cé 

1B D7 

1B D9 

BR FO 

8B FA 

BR CB 

Dl £0 

Di D2 

D1 D3 

D1 E0 

Di D2 

Dl D3 

D1 E0 

D1 D2 

Dl D3 

2B Cé 

1B D7 

1B D9 

D1 E2 

Di D3 

D1 EA 

03 C5 

63 D2 00 
81 FA 7FFF 
77 03 

EB 06 90 


40 
81 DA B000 


AS 0161 R 
89 16 0163 R 
D1 E0 

Dt 2 

oe 0? 

D1 E9 

D1 D2 

72 04 

FE C1 

EB Fé 

D1 EA 

D1 DB 

BS BO 

2h D9 


’ 
SUBTRM: 


UNIFOL: 


FNDEXP: 


EXPONT: 


e 
’ 


RCL 
RCL 
SUB 
SBB 
SBE 
FOV 
HOV 
OV 
SAL 
RCL 


RCL 


SAL 
RCL 
RCL 
SAL 
RCL 
RCL 
SUB 
SER 
SBB 
SAL 
RCL 
SHR 
ADD 
ADC 
CHP 
JA 

JHP 


INC 
BB 


MOV 
MOV 
SAL 
RCL 
XOR 
SAL 
RCL 
JC 

INC 
INP 
SHR 
RCR 
OV 
SUR 


Dy, 1 
BY, 1 
AX,SI 
DX, DI 
BX, CX 
S1, AX 
DI, DX 
CX, BX 
AY, 1 
DX, 1 
BY, 1 
AX, 1 
DX, 1 
BX, 1! 
AX, ! 
DX, 1 
BX, 1 
AX,SI 
DY, DI 
BX ,CX 
DX, 1 
BX, 1 
DX, 1 
AX, BX 
Dx,0 


DX, 7FFFH 


SUBTRM 
UNIFO1 


AX 


DX, BOOOH 


LOSEED, AX 
HISEED, DY 


AX, 1 
DX, 1! 
Beet 
AX, 1 
DX! 
EXPONT 
CL 
FNDEXP 
BX, 1 
AX, 1 
BL, BOH 
BL, CL 


‘BX, DX,AX, now contain Axx 


eget K=AIX/2°31 to BX 
sand X’=A8X(mod 2°31) to DX, AX 
scalculate X’’ = X’+K 


ssee if X’* ¢ # 
‘done if so Since X’’ can’t =f 


;else subtract N by adding 1 
sand subtracting 2*31 


sstore new x 


;Divide X by 2°31 to aake 

sa unifora {0,1). 

‘Set a counter to zero. 
‘Shift left and count the 
snuaber of leading zeros. 
sLeading 1 found. 

‘Not found, add { to counter 
sand look at next bit. 
sShift back to correct value 
sand suppress leading lL. 
‘Compute normalized exponent. 


; BL, DX, and AX now contain a single precision, uniform 
; (0,1) random number in BASIC’s floating point format. 
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O14E IF POP ps ‘restore segment register 


O14F BR 7E 96 MOV sD, CBP]+& = sget addr of variable U 
0152 8B 25 MOV (D1), AH ‘pass LSB of mantissa 
0154 88 35 91 HOV (DI]+1,DL  spass MSB of mantissa 
0157 88 75 02 MOV (DIJ+2,DH spass HSB of eantissa 
O15A 88 SD 03 MOV = (DI)+3,BL pass exponent 
15D 3D POP BP srestore BP for BASIC 
Q15E CA 0002 Relea sFAR return to BASIC 
) 

0161 4130 LOSEED = OW 4{30H sstorage tor 
0163 OQOAB HISEED =. OW OABH 
0165 RNGEN ENDP 
0165 CSE5 = ENDS 

END 
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Ber END PASE 


SCRNSHFT.SRT SOURCE-CODE LISTING 


cHEEESEEESESEPEEEEEEEEEEELEEEEEEEESEREEEEE EE EE 
*SCRNSHFT.SRT, Ver 2 -- BASIC USR Routine 
sby R.E. Nelsen -~ 27 duly 1984 


3 

sThis subroutine shifts 5 coluans to the left 
‘the 4 rows of pixels corresponding to the 
saodel's state vs tie graph. It is called 
sfrom BASIC using the comsands: 

5 

SCRNSHFT = 0 

DEF SEG = &HIA00 

CALL SCRNSHFT 


=e ‘we 


“2 =. 


‘This subroutines calls the ROM BIOS Type 10H 
sinterrupt (Video 1/0) using routines 12 & 13. 
cEEEEPESERELIEEEEEEEEPEEEELSEEESEEELEEEEEEE EEE 
; 
0000 CSEG SEGMENT 
ASSUME CS:CSEG, DS:NOTHING 


j 
sEstablish a file header for use by the 
sBASIC BLOAD comaand. . 


j 
0000 HEADER: 
0000 FD DBS s«OFDH sCode for BLGAD file 
0001 0000 \) *Seq addr location 
0003 0000 DW 0 ‘Offset location 
0005 0086 DW RIN_LEN ;Routine length 

§ 
0007 SCRN SHFT PROC FAR 
0007 35 PUSH BF sSave for BASIC 
0008 FB Sit ‘Enable interrupts 
0009 BY 0024 MOV = CX, 36 *Scrn coluan 36 
000C BA 0018 AGAIN: MOY «DX, 27 sScrn row 27 
OOOF B84 OD HOV = AH, 13 sRead pix at row DX 
0011 CD 10 INT {10H s& col CX. Pix in AL 
0013 83 E9 05 SUB’ Ek5 *Shift left 5 col’s 
0016 B4 OC MOV = AH, 12 ‘Set = 12 to write 
0018 CD 10 INT 10H spixel to new pos n 
OO1A 83 Cl 05 ADD CX,5 Ponift right a, col s 
001D 83 C2 10 ADD ODX,t6 sRom 43 
0020 B4 OD MOY = AH, 13 sset = 13 to read 
0022 CD 10 INT 10H sRead pixel 
0024 83 E9 05 SUB CX,5 SSiittelent eo 
0027 B4 OC MOV = AH,I2 « 3Set = 12 
0029 CD 10 INT 10H ‘Write pixel 
0028 83 Cl 05 ADD CX,5 ‘Shift right 5 col’s 
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002E 
0031 
0033 
0035 
0038 
003A 
O03SC 
O03F 
0042 
0044 
0046 
0049 
0048 
004D 
0030 
0054 


0036 
0057 
0059 
0058 
0050 
0060 
0043 
0069 
0067 
0069 
00C 
OOE 
0070 
0072 
0075 
0077 
0079 
0078 
OO7E 
0080 
0082 
0084 
0085 
0089 


83 
B4 
CD 
83 
B4 
CD 
83 
83 
84 
CD 
83 
B4 
CD 
83 
81 
72 


4} 
B4 
CD 
GF 
Bg 
8A 
B4 
BA 
CD 
83 
B4 
BA 
CD 
83 
B4 
BA 
CD 
83 
B4 
BA 
CD 
4} 
81 
72 


C2 10 
00 

10 

E9 05 
OC 

10 

C1 05 
C2 10 
OD 

10 

E9 09 
OC 

10 

C1 06 
FO OLA 
Bb 


OD 


Dg 
O15 
0018 
OC 

¥ 

10 

(2 10 
oc 
GS 

10 
C2 10 
0C 
C3 

10 

C2 10 
OC 
CS 

10 


F9 O11A 
DS 


;Now repeat procedure for regaining two rows. 


} 

ADD 
MOY 
INT 
SUR 
MOV 
INT 
ADD 
ADD 
MOV 
[NT 
SU8 
MOV 
INT 
ADD 
CMP 
J8 


’ 


Dx,16 
AH, 13 
10H 
Cx,5 
AH, 12 
10H 
Cx,5 
DX,16 
AH 13 
10H 
Cx,5 
AH, 12 
10H 
ers 
CX, 282 
AGAIN 


‘Row 59 


‘Row 75 


‘Shift to new coluan 
‘Is last col reached? 
sift no, do next col. 


sFinished with rows, so now aust blank out last 
su Coluans of each row. 


§ 

INC 
MOV 
INT 
MOV 
MOV 
NEW COL: MOV 
- MOV 
MOV 
INT 
ADD 
MOV 
MOV 
INT 
ADD 
MOV 
MOV 
INT 
ADD 
MOV 
MOV 
INT 
INC 
CMP 

J 


’ 


cx 

AH, 13 
10H 
BL, AL 
BY 277 
DY ,27 
AH, 12 
AL, BL 
10H 
DY, 14 
AH, 12 
AL ,8L 
10H 
DY, 16 
AH, 12 
AL, BL 
10H 

DY 14 
AH, 12 
AL, BL 
10H 

CY 

CX, 282 
NEW COL 


sCOle2b5 

sSet to read dot 
sRead blank dot 
:Save color code 
‘Col 277 

*Row 27 

‘Set to write dot 
;aet backornd color 
sWrite blank dot 
;Do again for 
srow 43, 


sand row 399, 


sand row 79. 


sNew col on right 
sLast coluan? 
;No, so do next col. 


;Done, so recover BP and return to 8ASIC. 
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008H SD 
G08C CE 
0080 

= 0086 
008) 1A 
O08E 


SCRN_SHFT 
RIN LEN 


CSE6 


POP BP 

RET 

ENDP 

EQU $ - SCRN SHFT 

DB O1AH ;Need for BLOAD 
ENDS 

END HEADER 


85 


BIBLIOGRAPHY 


Cheney, W. and Kincaid, D., Numerical Mathematics and 
Computing, Brook/Cole, 1980. 


Coffron, J.W., Programming the 8086/8088, Sybex, 1983. 


Davison, R.J., Graphic Simulation of the Poisson Process, 
Master's Thesis, Naval Postgraduate School, October 1982. 


Howson, H.R., "“POKEing Around in? Ehe LBM Pe, Pare 
Accessing System and Hardware Facilities," BYTE, v. 8, 
November 1983. 


Howson, H.R., "“POKEing Around in the IBM PC, Part 2: 
Developing Subroutines for BIOS Interface and Screen-Display 
Disk Storage," BYTE, v. 8, December 1983. 


IBM Corporation Personal Computer Library, BASIC, May 1982. 


IBM Corporation Personal Computer Library, Macro-Assembler, 
May 1982. 


IBM Corporation Personal Computer Library, PCDOS 2.10, 
September 1983. 


IBM Corporation Personal Computer Library, Technical 
Reference, July 1982. 


Law, A.M. and Kelton, W.D., Simulation Modeling and 
Analysis; MeGraw-Hili, igs2. 


Naval Postgraduate School Report NPS55LW/3061A, Random 
Number Package LLRANDOM, by G.P. Learmonth and P.A.W. Lewis 
June 1973. 


Rollins, D., "The 8088 Connection: Interfacing IBM PC BASme 
to Machine-Language Programs," BYTE, v. 8, July 1983: 


Ross, S.M., Introduction to Probability Models) Zameaua 


Academic Press, 1980. 


Scanlon, L.J., IBM PC Assembly Language: A Guide for 
Pregrammers, Bradyp i 7oee 


Stein, M.L. and Munro, W.D., Introductions to Machen 
Arithmetic, Addison-Wesley, 1971. 


86 


itv Droit hewn LIST 


No. Copies 


Defense Technical Information Center 2 
Cameron Station 
Alexandria, Virginia 22314 


Library, Code 0142 2 
Naval Postgraduate School 
Monterey, California 93943 


Colonel F.R. Giordano, U.S. Army 1 
United States Military Academy, Code MADN-A 
West Point, New York 10996 


Major R.E. Nelsen, U.S. Marine Corps 2 
Headquarters Battalion 

Headquarters, U.S. Marine Corps 

Henderson Hall 

Mealingeem, Virginia 22214 


Doceon 1.0, Schonreld 1 
Traffic Theory Department 

Bell Laboratories, Room 2M-626 

Holmdel, New Jersey 07733 


Premed. Di samy Code S5Ey 2 
Dept. of Operations Research 

Naval Postgraduate School 

Monterey, California 93943 


87 

















if eA Berney 
ie rhihe PaeA ERI oi ay, op Ria a ay My pga oe edie ao r i 
a 


ef Ub BEF, Tee Foti ws) ht a Fg nS 
oe 


am? ee Pee Fat 
mes at nny fey Be Cort &% ey ry if f M 4 
“aes of 


A656 od. 4" emia 3. an 4 
are ay ad Taber Aa LARD Ree 
Mi Tot4 “ eer 5 Y 
eens ott renter v ee ay ay 


“ ra Py rae 
pee tie te OO re oY ws 
; es Sear? Sar ET Pane iP | if pri 
arnt tT > Paty bt Meh 8 | +06 
oan fain Bae eis PP pt ait ed oy poy wy aA an Ce 
t 
i 







uf, 
A A 

trpphaeh it eo Mall ee aI C4 

SRR eA Ot Sit aa bara Ee ed 


ie ay OY ar ns 
8a mea iste’ 4 
oe a ree Se 
fi 


ee ote Oe 
ARB: Bm Pay FAC terete ta ey PAR) 
vm, ELT oe ae ear 


Fp Re pa Oe rea Srey 
Wed arta otdkat tan Aeed nite yen re | 
VF Prey wrong jetta dy yd Taree) led ee sehen re NG ft Sie Fla ap “f aD eestor <1 ee p wT 
OPS Mehrotra ort Sen Ee sw: ong hay igh oe ei 

7 ei etal AE ey Pre bd aM Rel ptr} | ET TN Ee Os Fagg 
HAE) ry ae ia pA ar eh as wo, Pe Fe eo pete de Per ty O75 eA te 
Pasith oa! rere ah Pear we POL at ane eye 




























Res oe FO ae a ee es 
Perea 1) iter ed PAG Vine hipal iff A, pees ’ ‘ een y 
he ce yey peters beeoded natal desi tegiy fro cred Ape at cr Ee 2am yy Pee Pete Pt ra rich ha Ab EX. 
y oath ie i Tatra PRENE G +: ao tty ae Ph itd MA AIA DMP SLAM the g Bas 4 
= eke oh Pn ks Fagin! rz Ce tM Tay ae pet Nd Kg $s Ps D 
tee ae POE Sadie pal no PET a Bee Aare Cs ie Pee! IF p a ee PN , eH Foy a LAS a 
Need Brae ft ‘ Sree lage. ry Pith seats ee STEN aT A OT Center *818 EGIL Emi ats e's 
Leite ees et eA Pep La $ es! Lie CeO ki LALA ahs HB Pee Tee ore Cr ye a Te i ee 
eres Cas peer Lye) me Pee edle it By Ks > yd at Pe LUE Lu FRA pre LASS Crane ne ee cy ew 
Perit Ain Phe ie wd Porn ee ine nets eee Phd om es See . Cin i HOT A 2 2) OR ree 
cod pre, are Prank x oe bg Febery yar we hte Pt “ MALAI 2 iF Sr ee is A Led Yt tl a 6 a ae 
5 ire rie he a erent i or Ce erp Py te ase fo taf yrs pratt jo8 putin a. i age Ae Ane eA 
Mikel mest Eee Ho Fitighss ibd EEO PP rewenntl anheded, Racha Rec a ke Dafa te ae cn ema 
oh VA per yr oe age nde rr el gg Speer A rd. Pea page cdl te hd vocal sel Seat oh La a I Par is 
epee a; ee Se a es Pe 2 byt SF, barat oe, ee ape RRR eo T a Se AS ee 
aay Dera Monhes lags hy hy hy harsh hy wae Ri PP ats ee ee RAAT PART aah , Es ee eT | 
eer ba ae rey i Ln yh Pa noah Ce an he hw atte ah ee Cy SURE on 3 
D Ad Bid Poo solar Nai dl be te Ce er ron Re aR et We tae he ashes ste aie , 3 n 
sa Pepa tee ae ce LS ares OTE OS ¥ aot # KP s Pe Ries ey “ee ao oy) eee tal OF oF Te ete ae 
Snr sees 2 2 fe ae. car eat ea One ore , Fiat Td pee Lyre Rare ie Mee ekg ) 7 
a Lh! ae Fw ae rt bh Ee eee’ ty Cy ae wee | Pararey its et Seka s a rae re 
mabe; FS Pee Ril Ca me ee | rien re ne aoe) Sg 14-6 ts p f. 3 < 
MEF od, 6 > ~ Th ee Be BY Sh L8 Beall toate cd vor ey a ioe | s) 4 
rea ap line he hia Tap P Ere E eo Seo py py ae Pe a © ota hte a 
is = 4) Asin $48 og 


huh ft ee erway a Cee oe eee Be ie 

2 eee CTU ay te a a 7 fs By! aT ae vi) Pi Fs 4 y LM 2 
AY y Apres Poy EE arse rb ae OTS ae Pee ay preemie ed ee yest A rar els Bei dea atte Fe de Lh | ote area Ae 
J J ete pr pe en eee * ig *y Me ts a tree ers covey i Ce Teer ars i PEN a Kate fh 

Tray ears ta38 fae 3 PED A® Rea ert ae Wwe had erate RCD Rosa Oe PR oP ee rae 
ted ay rae eeyy Pre pre, Apa bemin wrvicnitety Raat ne PO IP COL TR Oe he ee RK Bel ea ee LP ame iA mr) ae a 

eae Res H Ly PTE rtrd Ag auee : Y 4 
at rit ape ht pe ne ee me 3 oe Cor ws wrmrey: Re ee See s Seer AR ED et AE ry pay he raed ee eGo 
: PAO Seat EE Ue Paiyy A ooh | ad Tos eee ONY ee ay ies ee 

Pe ate Pee SO a a Bey Nk Rotel Spe roy. 





a ee ee Ses ea to , 














Rates diese fan OP hon ao Re Ar ree els Lamina! aa Me Ps 4 oes aE, wire: 5) i 4 fs rr ma 
re 3 J - Piette dot hives ee aed CO ent) oe ean eon Pr 2 he eae Peet rae ae ya 2 A 
ee peer peetha hee Pp ings PNT ep ee Se Nal a a2 ht ad ‘ teh a in oe ea ee r ape ¥; Pid ha ens ee & ied = 
fess mid gy Pre Foe OE T oe RRA AEP PUY , aE ew Oe Ye Bet vo ey Meany oo eae SO Cw ar ee 
f Apt rere: Latte ret ke a Care iY aw Yror es ee Se Pot 7a y 
3 Lat Per ae x CEL aera ars TL asap + 
ve & ean SNe sel 





Sir auy ¢ 
eat bor ie ee rE doy er a SCR Vee a “ie 
PUP Ca ae at a cae | dhe eg at LE Ee Ce ee ee a rt re ae 
wile es i dsi os Mita ld et th eek! Be PR se: Che oP Co kev Oop dw £4 ie a eee bc a 
oe vi phd a ocd Zeb beech tt OG ST oe ee ae ae | Alek 4.085 eG ea Sar eT ae Bee Mert eh a oe 
eM ie pects Pi rome one wlGh Cee Ll ehh oe ae ee re ary v6 be Rare ry iF eek Pm Oa ede | 
roa war) rf sate Fe balk alate eC Re Sere? ws CY arty 3 ean? ra en Ac Kaan A 
es wy é oh PP one indie ty Cita HP graben Gale h 4: thas Bion Sade Ped tes Sor Pad ety Vt Panera Sr hes Ags 
tao bien be sabimgant ies Fp feitecio air ylee Le ier a a ear abot a: NS es Seta boar SORE 
ppt ites Reb 4 add 4 pare netepehn nS PEPE EN TH a a 
Ba oa ie tk Fa ap ee Suna he as TE at Re el 























pate A Cada ee 
Phare hee ree TTT ieee Pe Wb ots A Baas PPS oes Pay ears wf : 
vee oh eae ERT Cop reR rere JB ce, > rb eee Agia J 4 rw E pghtte i) rr es a4 Va ie Pee ze rr Les BS 
oP eicy e he aaa Soh By a Cin coe) se ar Pee La ee eee ee ee rare « mr Par ye te 
ala el eo ra Cre en a cee 6 oer Fat Pie anh ale had Lol we Eto pe ee, eOee a pared aL E Oe haa ert ew Oe SR ‘ay 
See: amd et we et Sane Pint “ey i arp teh tod OR MELOY id SY re CE Pee ar ss ae Le eer) 
eae ae aren Bo eee Sis ee oe A bit dg beara UC Ie ee Pi BA 














PA eat 
eae rar s opt 









La Pek Ce te 
et: ele ee ee ae eee 





re 2 


OPA OG ed whieh ene a: ad 
yea le hes eee abet qi a id kN 





















hae ae | PeUueusL 
vs Pre Tae) Pacht Sale Lae pe Fa } Deen ee a eens at ely a 
bf ie ra ee i kta py oP ir. Cre KF opts Pre FF Li ee ar De be sic Pa Rs a z ore aS ft 
Fda dre : bese F x eer en a iva ik ieo5 pk See ee: OH agte de 

oe Fis | set: ie ere | Vik A ie fe See Te a eS Parner or gh, herr ae 

ra > bed ee rr rare ety Par a apn A oN mu oe Ng Parry, ete Kn eve 

IES oh pppoe Hel Z PET pe ‘4 ee Wile rN ere pee > A 
Pes rae Ts ee } CA Re ee 
A ae A aes ih phe: 


ey Cae ae) 
soiled gO ee este rw Ay sis it ped $ 
s “A MM ify pee me he sist weet ag fy hi 
We ere hes Parte the eer ee ; WP Te a pa 
gee Ei tiotat eee ia a ators se hee Or Peay uA td ba 
td ce ane | ig ‘) al 

Coed hoe Tet ry sia Sits Foe; ae 
at ee 2 Py a a 6%, 
ery 8 ke oe fe vy Le epee ee we 





har ees ey 

"ghee Aer. de gts g 5 tay a! at, Ages 
fe Wigtar, Pry ts bis\a-i 

ae by Se et ee | ; + 









try ae CP aes E = Tap ee os ee Oe eee fd wit AE dk OM Sg 5 
tae AN OER Oh ry datoleds i ope a ae any ie od aoe a rae ey 
Sal ad hat Ot aura ee pate Sale od, tasty saigid rae Fico! ca bea Rp noe Ae a ee Poi Pr vy 
Fob.b 8 Fock Deira) et ip et PING A LATS Te Ce ied Fo The Cr oT ae re yf OE ert i an 5 
Perey ere NPA Lae ea a TN SP ROR ere mae ETE ry er 
Pare Lr Labedah ase srt ad a ed ELE ers fj A , Xe 
ops a es Pree yoy a: Oi prerew Pe ha POE RA SEH i a “ea oe ip ce é ay, Pt ie een AZ mG r : 
ply et npr abst a Ra TET et oe EE eo 


hes sive po 
Be it rs TR AEE REST rari ¥ 
7 
tint phi 6 Rp. ery noe om rr 











i eter es eb ee ig 
a a 




















































s 
ar oY a te a fi) Py igh ms nl ee ate iat a 
A "Niel org Sem eee oe Ds Rhode «i Oe. be EE nt AO Wee tine a ae 
“ Cea re othe £4 teem aber ie eo ee a Sets Mee LE She a On hee : Th ae 
EEF AS SSA Syn AEE ES FSS Sed Pele Seite ene Tey aCe if se bare jected Ear ay fg are oP id Wik pobg et de Ge 
ries BET rates eRe Pt eres ap Leap cate ee iat Crue a ’ 4 a OP Ta BAS rs Pi eb k tne) Co os a oe 
E20 py He ak Ratan ap eae Fetes cS ody as ar $ee cy ee a oh ‘ale Ne La ord eae 
s ror H ee fs it de Ce gett injaake a or i 4 a) ota Yared Pe bw.’ peel bass 4, : AS carats f Pra ~ f fas 
fs eran aret SOE A KE tosh os eb. ct Ce Cee et PR: ole g thon ye et ae Tarr cf 
- pet h Saher Tike ee eae Oa ae ‘ et eee re Ff Ie ¢ 
aes ay rae Ae beta de nD al eM es Pr mea ae 
. s ess Sit: KF SEH § Crates Si ees ve ee Seta yi maar | 
ee Pipa oed Hears Pat} x eee Se Cre Re ee) fi A 
Pe a) Ste yee ee Eh Ets Re = 








We POT heat ae Pe ret 
We eres. ve eae ERD 7 
CES 3 tes a Pee ere boa ree ; 
oF ha ee ae) Lt ee Lat tat 
7 Pet eS Ue 

























Ay 
sees fi: 
aide ke 
She ta 










rc far Lo eee ae we an, 
Ie el Td is 


Pd Pe Ler 
a BTS ie tianiitc 
f UF a a ek 
Pare ay Stud rt ee Rae pre + P aa \ 

Fe i Re rasa Ca ae Yr a ie ; a fs gS Fae rats 
way iz es Tht ee Wer Pe at if ir ie a : re ot 2p 
Pai Bt z Pate we ot u rah F 1 of 
rr. SP aS a oF 





















Ts. 


“=e ve 


mr 


zy 









7 Leo. S 
Sunes Spud 5 
iG ts be 






























rs 
¢ A . * Pa 
pa es 1 a ‘ P a c : er 
yh eat 4 i 7 - Paes 
ered ee a ; ; . ae 
a tbe Silke? ¥ 
a a ~~ 
sisal fer 







: b, 
eh Pacis 33 adie wk a] Phiri ue 
ate LAS ate WS FRM eb 


Pee > “ 
Rika teetiel 
* i WAS ek iad CaS, eat iat 


A be 
eran eet te 
PRS Pre, A acetals’ Egle aa 


SRUR 















































































Le 
SP RR ee ea ok Pee x y yt eh : 5 Le 
a earch. es Bis we IAN ary curate RTineart HP Set tek COR eit thats ; : s 
re Wap Phe Se ka Zant fonlusstes te Md bere ee p Oh re rn pte bine Se a ees rare he : ei x Pe te ‘ T 
at og pret Sere ern tetae we MERE LPR ay ne, Sree t Sry ae ; WS Pe 
ph 01g ay phonation Soe a eae fipletd a a Bane. FT tore ae d ; tS 
mtagt hm ah ete is poate Ye ete hae ite ey >) 8 ay bir lhete 6 anit ah Be F or 
Sana Scere TA Ah bain ate aro: a oo a ; 
; “eye . an ui Hy a a fn we F 4 ! 
eee ae 4 ar ta SS AN SETS aa 5 igh on re ae ra Pa 
Se i by f 7% a hed a rv eas pia lle a me th tha 
St ee Ln A! ar Sr Tipsy Sark a ok da nip) At an Oy kad a repe x aes, " P E' 
Sle Rai GAs aaa ar kts ae + Fed Pats baat tetas 15 i aS var! eee ore +e I eee ae 
Li 5 Sty pa ad 4 Sei ery re eh aa ENE ah a7 re ¥. ae F ; F) tae he 
Hat pet De DN or ahi «Sag van Teo way >. orn ; ork 
4 WP * dye, MR ae & a's Q ASE ae .Y eal he Sone ye "4 WLS 
tae mad" eae fh eh he I ba PEERS Eat 2 | an! a me eh x ae ; s 5 ea 
thy = Deve arb! ee Sy tr perky S tea eo A ae ae Sy et eee ae ac ee A + SBA F 
atten th ALS lacs tt Wryens %\s or ut OS ee ~Om \ fe WW bk wy in YF 2 ee Pret Cb ee 
_ ae ars bi hae tate! SS" TOW Sten. Ere ey ry: ore Sins ha a fe - ae Ae etre Pee ‘ ts 7 
* a ay Boe ee See] ayn edie pe peat eee ae rete ee ta Pe cy a gd ae : ee shay a’, 
Weitere tn! Lert pigs Pe ss rir Soh oles xu) a Wheyg Wake aryeetisy ©. y : rh hk 
SN tte oar ts he te De ae pple » Ah wy ‘ y mir ae os SA Fe an 
Kae Pate: Sais “Es Eee ne at adeeb ony Se ee ee Ay <th pasveta es 
Sails gen Gee ett) TEA CENA oeghy “Hier, oss eRe oe 
m f hh EMyrheve Be 8 a c re ie h 
rent re hcntee * Pek panties Pere ah a RES Ry,! REN elie 
Sea Ore rath bbs peat bay sowta: ate) $s NE» iY. oe Tae SOT Pi 
ah cote tk ee Crate Ades Nae beta ad Ps te A | ee py Fi Pil oars VS as Oth a Cit Pee | 
" ears Re Tae ti tah a NEST: LP Bs a ea Se ay OES it Sra 
Soro reo Oo Aa A et Sa aC ACES es ay he 19 MINE NA Se ae a ct ar Kv. btw Sys 
pili Med ele data’ Pos oa hel Per eres . . Leah ht ie wee © teh rar a3 A Ce a ae tae hh a 
ie Won ene bt al phase Metis Sodas oS) Sri tart & yd Cian Ne la pes me wy: we bad ag , cr i Soe, ¥ €' 3 ae i, Ps ky ae 
“I oT) he , er .<, hey 4 Cia be Str es 
aa tt nthe eM ca ohreearenthines gue ote 3 os te oe Slee! oa i Be ho a a Rie 5 
Palio Se . 4 4 ; a 
SER TO eg arth ita de Adee ie deh tt ee oa: rah: ‘ 
hh Pan 


A kak Oh ae Pk 
er wf 





sv uy 





“ 4 " has ta a 
Por rR ad Pa eats. find Selb La Ae es es a] 








































oF L,Y 
tol he me he Sos 
aE oly ~ Goh! ith ‘ LE mn wn 
a uy a % nase 3 rei chr thd ra 5G {re SOs fe E Re cat SNS id OK A y \. “if Apc a> ae ape a 
ae AN are a Sa AS Reset ia 2 ae ee i ny Vw EA ey th Pee eat ro ke 1 Ce 
b a Caw iy) 5 a , 7 Lr : Sr : 
“ oe an Be ATR NERS ee aE ET ECR NE REGUS ON Sora kt OR SSPE A 
? Wed Sh he Sn KT Pate ; Pat , : ere 
plaed, gla. Nt re SANs Ses pes OO irc eee Coord aed NSS PAS ete tos POUR WS Ra eB 
OY to ao! Paeb) Sad RT my ie VARA 6.45 eA RIL A etiy Lees a ER rae ge gl RS 
ets Saath vs RA tS A eal a poeta ed athty bite h note ey aN ee CM aired Te WAN @M 7 1oN¢ ye 
rr a oy = es yi 28, oe i 
"i Se Rae Nk IS AT hee titel aah weet A ate hin ane. Ae er FON NCCI See see ae He ae 
RRS seer paige te Eu indo Wid det eta Laie Leh Th RMS ORR ir ice COCR rah Gur are 
_ Pe hs ayAd at NARS I NEN ONS le aie aan yas sla a aoe ta tae ohh ‘wire hy oe ad DATA P LA oe Me Hy 
a + re ‘ aaa Noa 1. ee CRT ee he Be oy ae: al be a ye tl dy Ae Py er Brat oo ory Ce ee Le a rae A 
Ae hte oe tir eX RLS Soe Pi eg: A BENS eal AL Pa aA ae Areey ame ah Ah On De 1 At SS mi 
. , : : on) ea nr "ee or c.9t b ' meant 
We G8 vectra ee agen ee ahha cheba oe Seve SiC Ns aR ah PE tk Mei So 
eee Tet SS 3 rh Mab date tat Stee dW oaks Pia a ar tia ye kU Be Sty Pee SOUR ey: reel 
a tate tek ee Te Or, UW R Ratbone ara hpi La 2s SE SO APOE A th Ae ay roe PTR RC ear 
SORE Eh, a 9 ihe Pocteruners Prag aT, ya ESS SW le ‘e's a ek Re oe hee ne .* a ee 
ANSARI Tt te hh ap PTA Se! PIRES 1H by ark Anaahat OE Ree Rita 
fat poole tia es bead. One a) ea 








a L 
& RTT G arero eee PS 7h. he ee ae oe 
1 OE yeaa OO Sak ee aus aay Paar eke ee Oe ee ee y 44 
PN CON hehe de Mh tae eh OR eS : 
8 





let Roh tt | hi Weta thins ey “ih ety At D 
aes ey Panto i ont “Ng 







































haa Po ee 


a aye * RD Be She ee Perv ey rene re re Se ni ge 
at enue CR Or eee Ae ea PORE Ra 544 Sta waa, AA aU ROLE IC Rb 
ee SAE mt Sp pice ha aera Ryn Say Na ote Ae Be nt ANG A ah rete FES ESS Marte 88 weet 
Sasi SAUER SEY ea a yoy par fpd bb teianL ES A rns aa nt ri far Sah Bd RK Ci +e Me Ney on aN aes, 4 : wi a 
ue i No rey} ee ys Cae Way tl Ge ew Corts piers . a“ F re rk * 
2 rama tk GAA BST Sue oi a ae eee SRT , ty " pry oo mi ra id 
Bas ke eet bes ~, = H oy Ph Ah 0 TO A bk Hebd te Ena ua Mier er we ark ata 
ae re Pet aA SROs etree ea Hi ae a A cabs BI see ea a St Ys ee EL Sora rire ts a re LV 5 by, 
mh Ok Aird — A ln ee wh rls amiga tale CODE ark tae “* ey J fest ts PL Ps rae he oe cae im He We as i 
eg un. hates le Mla a sa Neln aii he) SA aN aot) hs ta tot gh Hi, 5 i GR Le a CP: ‘ h EPL het 
‘hee SAR Mite Tata aot Nessrereesg a Noa OO Seah ar ae COAT ACR. BACK ty atria Rat lh cel Siam ro ee, 
SRS: BeeKt SME MERON TL TRESS EE IE AA OSA a ho ja SRN CL 
SRNL Sk ry REAR ARR rae a AN at pays ae ny MyAOL Xf) AC: Da ae ak Or 
Cake ast nee area EE Arora dh CH RO aL EY Unb thd PO WeCROrr orton arias Vr ‘ ee 
Pad tS DOE heeb RL bets 2 oe are Ms yh Xo i Pht Pr Lye gy ee ath a 
a] Sea hes Paty blac A etd i dhe Lie Pht a aR ly vs SEERA pd ot De es “of a tw « iy 2 ai 2 
nary Perens se Lar me RNs ry iA: he arte oy a gS WT WR tba deo et Ur ee CL 
ka oa Seam SEVER A ny Cary ma 4 MARL XS pe ae As rts Loe a Nid SU Sr Pek eh be A Se ee a 
: hese 4 : 
Tr a RRNA ia ri eC RRS AT Wea R a 


+ hate Lt mt A at hl 










D Rewer ers s , 5 Fa 
SEL Pe get yt ey SE UK ye 







es A Se WASP a Li ‘ ? 
é \ Gale br Sh RS ee | et Aree oa Beat oh a a ey 3 bP eeyrm a2 Qa 
Abilene techs ON RR Te § ae Peer " hk AE Va ae Ve RPE «ord be ba 
rt Q ef Tai Ce ee ee eereee Tey 1 ew 
AD iat WRG ET. rte ae IR Er abe He Be 


Lr ae ie Oe es 





Ce + 448 ¢ 


VTERII 
ed Po Bt ie 

a ae 
, open, 
Saw ae 
Uae ae 











hf 7 
. ra 





SC ve Agta ae 







Graphic simulation of a ae repairma 
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